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Abstract 



CZ ^ We describe how to use zpl Tj^ji, a program based on a semi-analytical approach to fermion pair 
production in e^e~ annihilation and Bhabha scattering. A flexible treatment of complete 0{a) 
QED corrections, also including higher orders, allows for three calculational chains with different 
realistic sets of restrictions in the photon phase space, zpl Tj^ji consists of several branches with 
varying assumptions on the underlying hard scattering process. One includes complete 0{a) weak 
loop corrections with a resummation of leading higher-order terms. Alternatively, an ansatz inspired 
from S-matrix theory, or several model-independent effective Born cross sections may be convoluted. 
The program calculates cross sections, forward-backward asymmetries, and for r pair production also 
the final-state polarization. Various interfaces allow fits to be performed with different sets of free 
parameters. 
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1 Introduction 



There is a growing demand for flexible programs to fit the very precise data on fermion pair 
production from experiments at the e"'"e~storage ring LEP I: 

e+e--^//(n7), (1-1) 

including Bhabha scattering, 

e+e" — ^e+e"(n7). (1.2) 

It is important that such programs allow for model-independent and Standard Model [Q fits 
to the data. In addition, it is interesting to be able to fit the data to theories that go beyond 
the Standard Model. Because experimental cuts tend to be more complicated than can be 
realized with semi-analytic programs, typically Monte Carlo programs are used to correct for 
such cuts and detector inefficiencies before fitting. 

In this article, we describe the subroutine package zfI Tj^^i- This program is based 
on a semi-analytical approach to the radiative corrections that are needed for the analysis of 



reactions (|1.1|) and (IT^). 




Figure 1: Scattering angle d in e^e annihilation. 

T 

The zpl Tj^R package employs an approach which relies on formulae that are analytically 
integrated over a finite angular region with respect to the scattering angle, as shown in fig. [l|. 
The program directly calculates predictions for observable quantities and not corrections to 
Born approximations. The total cross section, ctt, and the forward-backward asymmetry, 
ApBi may be calculated in a non-symmetric angular interval, ci < cosd < c^'- 



crr(ci,C2 



dcos^ 
ci a cos v 

apBici, C2) 



ot(ci,C2) 



where 



(^Fb{Ci,C2) 



C2 



dcosi!} 



Cl 



da 



dcosd 



These expressions are constructed from the following integrals: 



1.3) 

;i-4) 

;i.5) 



a(0,c) 



a cos 17 



dcos^ =2[-Hc) + a,,(c)]. 



;i.6) 
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Here o"t(c) and <7fb{c) are two special cases of (Ja(ci,C2), which will be the basis of the 
discussion in the following chapters: 



c rO 



d cos 1!} — d cos 'd 



da 



dcos^ 

By simple algebraic combinations of the above constructs, one may derive various measur- 
able cross sections and asymmetries. One must, of course, take into account the possible 
beam polarizations and final-state helicities within the hard subprocess description {(t^^fb) ^ 
explained below. 

For reaction ( |1 . 1|) , excluding Bhabha scattering which will be discussed in section both 
functions ax and apB may be split into different contributions from initial-state radiation, a™', 
final-state radiation, a^'^, and their interference, a'"* {A = T,FB): 

aAic)=aT{c)+afic) + a\^ic). (1.9) 

Common soft photon exponentiation for initial- and final-state radiation, which relies on a 
more compact (but also more sophisticated) formula, has been realized in zpl Tj^f[- 

aA{c)=a't+''^{c)+o'l\c). (1.10) 

Alternatively, the program allows the user to choose a simplified treatment of the (small) 
contribution from final-state radiation: 

aA{c) = af{c) (1 + l^Q}) + af{c). (1.11) 

rri 

The expressions introduced in (|1.10|) and in (|1 . 1 1| ) are realized in zpl Tj^ji as one-dimensional 
numeric integrations over a photon phase space variable s' . 

Photonic corrections to the cross sections and asymmetries are implemented by convoluting 
the Born cross sections (cr^") with radiator functions (-R^): 



a\{c) = ^5Je /"^ dva^X^is, s')R'X{v, c), (1.12) 
dA Jo 



dA 

where a = ini,ini+fin,int; dx = |, dpB = 1; -s' = (1 — v)s] and v is the energy of the radiated 
photon in units of the beam energy. Further, cr^ contains the dynamics of the basic process 
to be studied, and the functions Ra depend on the treatment of the QED effects. There are 
several ways to describe zpl Tj^^. It contains: 

• three calculational chains with a different handling of QED corrections plus the Bhabha 
chain, 

• four branches which differ by the theoretical description of the hard scattering process, 

• seven interfaces with different choices of input/output parameters. 



5 



1.1 zfI T^r Chains 



No cuts - a fast option 

In this chain, the cross sections are calculated with formulae that assume that there are no 
cuts applied to the photon phase space. 

Phase-space cut on the minimum invariant mass of the // pair 

The underlying formulae may be found in l^,^ and in references quoted therein. This chain 
allows for a cut on the minimum invariant mass of the final-state // pair, which can be 
reinterpreted as a cut on the maximum of the allowed energy of the bremsstrahlung photon. 

Cuts on energies and acollinearity of final-state fermions 

This treatment of the photon phase space follows the basic lines of that of the above chain. 
The restriction on the maximal photon energy is replaced by a simultaneous cut on both the 
energies of the produced fermions and on their acollinearity P|. This chain also allows the 
calculation of differential and integrated cross sections for Bhabha scattering using the BHANG 
package which has been incorporated into zf^"^^^- 

Furthermore, in both the latter chains one can impose a restriction on the maximum 
production angle of the outgoing antifermionQ. 

T 

1.2 zpI T^ji Branches 



Analytic Standard Model formulae with higher-order corrections 

This is the central branch of the program. The calculations of the partial and total Z and 
W widths follow and [§] respectively. The explicit formulae for the improved Born cross 
sections with electroweak corrections are described in [^. In addition, improvements have been 
realized in the program by including various higher-order corrections, which will be described 
in detail later. The electroweak loop corrections are determined in zpI^Tj^^ using the DIZET 



package |T0[ for all channels including Bhabha scattering [11 



Model-independent ansatz using eff"ective couplings 

This approach assumes that the effective axial-vector and vector couplings of fermions to 
the Z are real, constant, process- and energy-independent as in ||12|, |13[. It is known from 
comparisons with Standard Model predictions that these assumptions allow for quite a good 
approximation. 

Model-independent ansatz using partial decay widths 

Following general arguments of field theory, one can describe resonance scattering with the 
help of the partial decay widths of the resonance. This is particularly advantageous since 
measuring the Z line shape allows for a very precise determination of the partial decay widths. 



A compact description of the underlying formalism may be found in [12 



S-matrix ansatz 

The cross section ansatz due to general S-matrix ideas as described in has been imple- 



mented in the program. The main advantage of this branch is that it gives the mass and total 

^ As a matter of convention this cut is imposed on the antifermion only. Because of CP invariance the cut 
could equally well be applied on the ferniion instead. 
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width of the Z with minimal assumptions on the underlying dynamics. As with the other 
branches some additional degrees of freedom are available; however, the physical information, 
which can be extracted from them is limited. 



The advantage of the various model-independent branches is the simple picture of the dy- 
namics and the gain in flexibility compared with the Standard Model. Such model-independent 
approaches also allow for pragmatic checks of Standard Model predictions and practical gains 
in computing time. If different branches give statistically significant differences with respect 
to the various parameters then perhaps a strong indication of New Physics exists! 

Thanks to the flexibility of the convolution approach to QED corrections in zfI Tj^^, it 
is relatively easy to make different assumptions on the hard scattering process. While the 
above branches cover some of the most important theoretical tools for LEP I physics, it is a 
straight-forward job to add new branches so that predictions for New Physics can be made 
within the zfI^Tj^^ framework. One such example is described in where the mixing of the 



Z with an additional heavy Z' is implemented. In addition, other possibilities, which cover 
some New Physics by extensions of the weak form factors will also be discussed. 

1.3 zpI T^ji Interfaces 
Subroutine ZUTHSM 

Calculation of Standard Model cross sections and forward-backward asymmetries as functions 
of Mz, rrit, M^, and a^. 

Subroutine ZUTPSM 

Calculation of Standard Model r polarization, Apd, and r polarization forward-backward asym- 
metry, A^j^, as functions of Mz,mt, M^, and a^- 

Subroutine ZUXSA 

Calculation of model-independent cross sections and asymmetries as functions of the normal- 
ization form factors (p), effective vector (v) and axial- vector (d) couplings, respectively. 

Subroutine ZUTAU 

Calculation of model-independent final-state polarization in r pair production as functions of 
the normalization form factors, effective vector and axial-vector couplings. 

Subroutine ZUXSA2 

Calculation of model-independent cross sections and asymmetries as functions of the squares 
of the normalization form factors, effective vector and axial-vector couplings. 

Subroutine ZUXSEC 

Calculation of model-independent cross sections as functions of the partial (P/) and total Z 
widths. 

Subroutine ZUSMAT 

Calculation of model-independent cross sections, based on an S-matrix inspired ansatz, as 
functions of Mz,Tz, etc. 
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The above interfaces have been designed with the analysis of LEP I data in mind. In fact, 
the accuracy of the Standard Model branch of zpl Tj^^ has been optimized near the Z pole. 
Nevertheless, the Standard Model branch of the package can be used at PETRA, TRISTAN, 
and linear collider energies without changes. Many examples of the use of the zpl Tj^j^ package 



exist in the literature (see e.g. |]T6|. We will thus make no attempt to describe how to use 
zpI^TjT^R to fit data. 

The organization of the article is as follows: section |^ describes the treatment of photonic 
corrections; the description of the various theoretical treatments of the hard scattering process 
is given in section |^ for the Standard Model branch and in section ^ for the other branches] 
the search for effects of New Physics with zpI^Tj^^ is commented in section ^; initialization is 
described in section 10; section documents the interface structure; finally section M compares 

rp 

zpl Tj^ji results with those of other programs for weak mixing angles fl^, |18[ and widths [|1^, 
|20[| , and also cross sections in an energy range covering both PETRA and LEP I energies and 



beyond (ZSHAPE |21, 22|, ALIBABA [^3]). The appendices contain a description of the contents 
of some of the common blocks of zpl Tj^^ and an example of the use of the package. 



2 Chains of zpl T^^: 

Photonic Corrections with Different Cuts 

In this section, we will describe the functions (Jt{c) and apBic), which were introduced in 
( |1.7|) and (|1.8| ). A complete treatment of photonic corrections would also include the running 
of the electromagnetic coupling constant, which will be discussed in section ^ on the hard 
subprocess description. 

In order to get a finite, gauge-invariant result, real photon bremsstrahlung from the dia- 
grams of fig. 1^ has to be combined with photonic vertex corrections of fig. ^ for initial- or 
final-state radiation and for their interference with the box-diagram corrections of fig. |. 




Figure 2: Real photon emission from initial (a) and final (b) states. 

T 

zpl Tj^j^ relies on the following numerical integration for the contributions introduced in 
( |1.10| ). The common soft photon exponentiation of initial- and final-state radiation is taken 
into account with: 

<'+''^(c) = ^3?e dva\{s')RX\v, c)Rf{v), (2.1) 



A 



where s' = (1 — v)s. Final-state radiation is described by R^, which is more complex than 
a simple angular integral of R^{v,c)] in it has been shown that R^ is almost completely 
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Figure 3: The photonic vertex corrections for the initial (a) and final (h) states. 



7 

Figure 4: Box diagrams with virtual photons, which combine with the initial-final interference 
bremsstrahlung into finite and gauge-invariant contributions. 

angle- independent. For each of the cross section parts, the contributions from 7 and Z ex- 
change and from their interference can be separated: 

^Ais) = E ^Ais; m, n) = 7, 7) + ^l{s; 7, Z) + Z, Z). (2.2) 

m,n 

For the interference of initial and final states, this decomposition is unavoidable. This is due 
to the differences in the 77 and 7Z boxes in fig. ^ which regularize the infared divergence: 

aj\c) = -^3fJe dv J2 ^^^(s, s'; m, n)E}2\v, c; m, n), (2.3) 

a A Jo rn,n 

where m,n = 'y, Z. The origin of the complex structure of the initial-final interference brems- 
strahlung contribution is two-fold. First, the cross section part originates from the interference 
of matrix elements with emission of a photon before and one after the hard-scattering process. 
This leads to the dependence of the hard-scattering cross section, a^, on both s and s'. Sec- 
ondly, the virtual corrections of initial- and final-state radiation or of the interference have 
different structure (fig. |^). The simple vertex diagrams (fig. §) of the former factorize into the 
Born cross section and a universal factor, while the box diagrams (fig. |^) with two-particle 
exchange do not. This leads to a dependence of the interference radiator functions, -R™*, on 
m, n. 

The radiator functions i?^(f,c) are the result of a three- fold analytic integration of the 
corresponding photon phase space: 

i?^(f , c, m, ra) = / dv2 / dcos-^} / x^(s, f , ^2, cosi?, <^-y), (2.4) 
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where Xa the result of a Feynman diagram calculation. Further, s' = Rs = (1 — v)s = Mj^ 
is the invariant mass of the fermion pair, V2 = M|^/s and (j)^ is one of the photon angles in the 

(7, /) rest system. Two treatments of the photon phase space are realized in zfJ^Tj^^. These 
are shown below in the Dalitz plots of figs. ^ and 0. The variable f 2 has been integrated over 
analytically, while R remains to be numerically integrated by zfI Tj^j^- Note that the corner 
of the photon phase space, which corresponds to the emission of a soft photon is located near 
R = l. 

As may be seen from the definitions ( pTT]) and (pTSp, the angular acceptance cut, ci < 
cos-i? < C2, limits the scattering angle "i? of the final-state antifermions (see fig. [l|). In this case, 
the scattering angle of the fermion / remains unrestricted if the other cut(s) do not imply an 
implicit restriction (see section |2.3|) . 



In zpl Tj^n^ the QED contributions include the complete 0{a) corrections and soft photon 
exponentiation. It should be mentioned that the radiator functions (flux factors), R% differ 
for different observables [A = T, FB) and for different bremsstrahlung origin (a = ini, fin, int). 
In addition, the radiator functions for the integrated cross sections and the differential cross 
sections are not the same 0. Only at LEP I, around the Z resonance, do all the radiator 
functions agree approximately p4|, Some other semi-analytic programs use equal radiator 
functions for the total cross section and for ApB- At LEP I energies, where hard photon 
emission is suppressed, and for loose cuts (thus not enhancing the initial-final interference 
terms), this is numerically acceptable. zpI^Tj^jii however, uses the correct radiator functions; 
the underlying formalism thus allows an application of zfI Tj^pi at energies far away from the 
Z peak. No part of our treatment of the bremsstrahlung is specific to physics near the Z 
resonance peak. 

Higher order QED corrections have been implemented in zfI Tj^^ for initial-state radiation 



contributions, besides the above-mentioned soft photon exponentiation, to as in |21] and 



to AiT^ as in |2^. For the two calculational chains, which involve an angular acceptance cut, 
these higher-order corrections are treated with an approximation that assumes a Born-like 
angular behavior. 

When no acceptance cut is applied, c = 1, the expressions for 0"t(1) and Afb{^) approach 
well-known formulae for ax H^, HTJ and ApB [il, • 



2.1 No Cuts 

rp 

Of the various calculational chains contained in the zfI Tj^^ package the simplest to describe 
is the one where no cuts are allowed. This chain has been realized with special formulae in 
order to make it as computationally fast as possible. Here the photon may have any energy, 
A, up to the kinematic limit: 

A = E^V^beam, A < A^^^ = 1 - Am)/s. (2.5) 

Thus, the radiative corrections depend on the fermion masses even for the light quarks and 
leptons. This dependence can be important when total cross sections are determined from 
experimental data, and is of special importance when comparing results from other semi- 



analytic programs. The latter will be discussed in section 
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2.2 Convolution Integral with Cut on the Invariant Mass of the 
Outgoing Fermion Pair (s^^j^) 

zpl Tj^ji allows for a constraint on the minimum allowed invariant mass of the outgoing fermion 
pair, s^j^: 

4in = (1 - A)s. (2.6) 

This is easily re-interpreted as a cut on the maximum allowed energy of the bremsstrahlung 
photon, A. In this calculational chain the s^j^ cut may be combined with an angular accep- 
tance cut: 

Ci<cosi9<C2. (2.7) 
The Dalitz plot shown in fig. ^ corresponds to a A cut of A = 1 — i?™™. 




Figure 5: Dalitz plot for the photon phase space with a cut on s' . 

The allowed region is a triangle in the ultra- relativistic limit. Note that f 2 is not influenced 
by the cut. This simplifies the analytical integration over V2. 

Explicit expressions for the radiator functions -R^(f,c) discussed above may be found in 
the literature. For initial-state radiation and initial-final interference, they may be found in 
eqs. (8), (18) in Q, respectively. For final state radiation, the angular dependence is relatively 
simple and eqs. (132-134) in |^ are valid. The radiator functions for common exponentiation of 
initial- and final-state soft-photon emission implemented in zpl Tj^^ are derived from eq. (157) 
in 0, as has been described in section 4 of 0]. 

2.3 Convolution Integral with Cuts on Fermion Energies and 
Acollinearity (^min^^max^ 

As an alternative to the s^j^ cut, one can apply another set of cuts on the outgoing // pair p. 
Cuts on the minimum energy, E'™''^, and the maximum acollinearity, ^^'^^^ of the // pair in 

addition to angular acceptance cuts have been implemented in zpl Tj^^. 

Figure |^ shows a Dalitz plot of the allowed phase space for the two energy variables R and 
f2, introduced above. 
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R^R, Re R 



Figure 6: Dalitz plot for the photon phase space with cuts on and ^ as explained in the 
text. 

The boundaries of the allowed phase-space region are defined by the following conditions: 

vT"- = 1-Re, (2.8) 
vf'iR) =Re-R, (2.9) 

^ ^""'^ - {1 - R,y + 4R,V2^ ^'-'^^ 



where 



2Ef- 2Ef- 1 - sin(r-/2) 

^«-l + sin(e-V2)- ^^-^^^ 

The absolute minimum of R is given byQ 



Further, the upper integration limit in ( |1.12| ) becomes 



A = l-i?™^ (2.13) 

The above relations are independent of the scattering angle and are, thus, compatible with 
an angular acceptance cut: 

ci<cost9<C2. (2.14) 
The turning point, Pt in fig. 0, of the acollinearity bound of the integration region is: 



Pt=[RuV2,t] 



R^', 2(1 ~ ^c) 



:2.i5) 



This is significant since it allows the user to apply a reasonable approximation of the acollinear- 
ity cut in terms of the simpler A cut; this can be achieved by using for the definition of 



the integration limit ( p. 61) in section ^ 



2sm(r"/2) 

* * 1 + sin(^'^^V2) 



^The current implementation of zpflTgj^ assumes that i?™™ i?™' 
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The quahty of such an approximation depends critically on the values of the -Emin cut and the 
^max ^^^.^ loose cuts it improvcs. 

Because of the approximations that have been implemented in the QED calculational chain, 
the user must be cautious in applying severe cuts. Since the approximation is ultra-relativistic 
one should restrict oneself to the region of the phase space: 

E"''">m/, (2.17) 



Near the turning point Pt introduced in ( |2.15| ) the validity of the soft photon exponentiation 



approximation comes into question^. To avoid any such problems the following restrictions 
should be observed: 

E°^^°<0.95Ebeam, > 2°. (2.18) 



This last limitation guarantees that the second-order terms [/51og(l — -Rcut)]^ with (3 = 2(a/7r) 
[log(s/mg) — 1] and _Rcut = Re, Rs, do not become too large: 

|;3l0g(l-i?eut)|<l. (2.19) 

This corresponds to 

« exp(/? p -l ^^^^^^ ^ expi-(3-'). (2.20) 

exp(/?-i) 

Finally, we would like to point out that the acollinearity cut has an indirect influence on 
the acceptance cut. It is easy to see that the maximal scattering angle of the second fermion 
(which is unrestricted by the user's acceptance cut) becomes limited by an acolhnearity cut, 
i.e. the scattering angle of the second fermion is limited to [— (^™'^'' + ^max-j^ j^^max _|_ ^max^j_ 

2.4 Photonic Corrections for Bhabha Scattering 



The Bhabha scattering cross section (|1.2|) arises from the sum of s- and t-channel exchange 
cross sections and from their interference. The s-channel part needs no further comment since 
it corresponds completely to ordinary fermion pair production. In the t channel the energy 
variables, which correspond to {s,s') are {t,t'), where: 



t = -is(l -COST?) (2.21) 
2 



and 



s' 



t' = t-. (2.22) 

s 

The t-channel propagator for a massless photon is proportional to 1/t or 1/t'; it thus becomes 
divergent in the forward direction, i.e. as — > 0: 

^^Bhabha I 

Such a divergence is common in calculations of the Bhabha scattering cross section and it 
prevents a reasonable definition of a total cross section without at least an acceptance cut, 
even at the level of the Born approximation. 

more advanced exponentiation procedure 129] circumvents these limitations. 
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For large-angle Bhabha scattering > 10°, this problem is absent. Near the Z peak, 

such a condition guarantees that the photonic t-channel exchange contributions, including the 
QED corrections (with an effective t' = s'{l — cos-i?) < t in the hard-scattering process), are 
at most of the same order as the non-resonating terms from the s channel. Of course, an 
acceptance cut does not prevent t' from becoming small because of the emission of a hard 
initial-state photon, in which case the t channel dominates. This divergence can, however, be 
circumvented by excluding very hard photons from the observed cross section with an ii^™™ 
cut. In any case, the hard photon corrections to these contributions must be carefully taken 
into account. 

At LEP I energies, terms with Z exchange in the t channel are strongly suppressed owing 
to the form of the Zpropagator (~ l/(t' + M|)) and contribute less than 1% to the cross 
section. In summary, the contributions, which arise from photon exchange in the t channel, 
compete with those of the s channel; however, near the Z resonance it is clear that the s 
channel must dominate. 

An explicit description of the QED corrections to Bhabha scattering which have been 
implemented in zpI Tj^^i will be presented elsewhere In order to discuss some features 
of the present implementation, we give an explicit example for the general structure of the 
Bhabha cross section: 



d cos d d cos d 



a Viy2 



Here the first term corresponds to the s-channel part. The sum under the integral extends 
over a, denoting in the s-channel diagrams initial- and final-state radiation, in the t-channel 
up and down radiation, and in the interferences the corresponding combinations. Further, a 
sum extends over (Vl,V2), the possible combinations of propagators 7^,74,^5,^4 from the t 
channel and the interference. In ( p. 241) , all functions i?" have the form 



R''[v,cos'd]Vi,V2) = S{1 
where 



(X 

l + -S^{s, cos^9;\/i,\/2) 
7r 



Af^ + -H''{v,cos^;Vi,V2), (2.25) 

TT 



1 /3 = 4^ log^ + losi^ , (2.26) 

TT V mi 1 + cos V I 



and i?™™ was introduced in ( p.l3[ ) 



The functions 5" in (|2.25| ) contain the soft photon (plus corresponding virtual) corrections, 
and if" the complete 0{a) hard photonic corrections. In the s channel, the hard photon part 
depended on s'/s only, while here, due to the t-channel propagators, it is also dependent 
on t'/s or on t/s. As a consequence, it looses its universality and depends also on the kind 
of bosons which are exchanged ( i.e. 7 or Z). Further, the running QED coupling (if not 
assigned formally to the radiator functions, it is contained in the hard-scattering cross section 
0""'°) depends, in the t channel, on the scattering angle as well. 

These {t, t') dependences have the far-reaching consequence that the integrand in (|2.24|) 



depends in a more non-trivial way on the scattering angle compared with the s-channel case 
- thus preventing an analytic integration over cos -d, which is the basis of the fast computing 
of zfI\r. 
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In the current implementation of BHANG, the cut conditions of section |2.3| are taken into 
account in the functions ( |2.25| )f|. 

Further, in (|2.25|) the leading higher-order corrections due to soft and hard coUinear photon 
radiation with t-channel participation are taken into account in an approximate way. The 
cross section of the hard process is considered to be independent of the actual energy scale, 
i.e. assuming s' = s and t' = t. At LEP I, the error induced by this is definitely less 
than 1%. A simple improvement could be the choice of some better scales for the effective s- 
and t-invariants in the hard cross section, which effectively take into account the change of 
kinematics due to radiation. 

Hard photon radiation is considered in the collinear approximation for the cross section 
parts which correspond to Z exchange in the t channel, i.e. the appropriate functions 
H\{y, cosi?, Zt, Zt) are set to zero. 

The user of zpl Tj^ji should be aware that the Bhabha cross section returned with the aid 
of BHANG is to a much larger extent adapted for LEP I physics than for the other fermion 
channels and contains more approximations in the treatment of the QED corrections. 

part. 

3 The Hard Scattering Process: 
(I) The Standard Model Branch 

We now describe some general features of the cross section formulae for the hard-scattering 
subprocesses. In all branches of zfI Tj^^, we can denote: 

a^(s,s';m,n) = J^lm, n; s, s')^ [/C^(s')/C* (s) + /Cm(s)/C* (s')] . (3.1) 
For initial-state radiation this simplifies to: 

(T^(s', s'; m, n) (y^is'] m, n) = Xyi(m, n; s')^m(s')^* (s')- (3.2) 
For final-state radiation s' has to be replaced by s. The propagator functions )Cn{s) are: 

S - Mi + iMnGn 

Here, A1„ are the masses and Qn are the widths of the intermediate gauge bosons. 



In addition to the QED-corrected cross sections ( |1.12| ), zfI Tj^^ can also return (improved 
or effective) Born cross sections, 0"^°™. These are constructed out of the expressions introduced 
above: 

aT^'is, c) = Da{c) {Ta(7, 7; s) + Wie [1a{i, Z- s)/C^(s)] + Xa{Z, Z- s) \1Cz{s)\^} , (3.4) 

DJA-{ 2(c+|c3) fora = T 

^A(cj- {oiA = FB. ^^-^^ 



* At present there is a limitation on the allowed value of the scattering angle i?; it must be larger than the 
acoUinearity This is due to purely technical reasons and this restriction will be removed in successive 

versions of the code. 
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The functions Xa contain the underlying dynamics of the hard-scattering process. Often, but 
not necessarily, they are assumed to be inversely proportional to s, s' . The different branches 
of zpl Tj^ji rely on various assumptions with respect to X4, as will be discussed later. 

For the photon (n = 0) the propagator becomes /C^ = 1, while for the Z, various possibili- 

ties exist in zpl Tj^ji- In recent years, much influenced by the discussions of the 1989 workshop 
on physics at LEP I |^|, it became common to use the following definitions: 



^Z{S) = ■ ,n , ... ^ , (3.6) 

s - Ml + tMzGz 

Mz = Mz, (3.7) 



r(.) - (3.^ 



where Mz and Tz are considered to be the mass and total width of the Z. This point of view 
reflects the fact that in a quantum field theory such as the Standard Model the Z width is 
predicted as a result of quantum corrections (self-energy insertions) and is, thus, naturally 
s-dependent. This s-dependence of Q becomes important only because the very narrow Z 
peak may be scanned with extreme precision, leading to errors of a few MeV for mass and 
width of the Z. 

The definitions ( |3.6| )-( p^ ) may be related to an alternate resonance description, which 
assumes a constant width: 

^^(") = M^'-M r ■ ^3.9) 

s — M| -I- zMzTz 



The following equality holds as long as the approximate relation in ( |3.8|) may be considered 
to be exactg l^: 

G^ICz{s)=G,ICzis), (3.10) 

Compared to (|3.6|) - (|3.8| ), (|3.9|) corresponds to another ansatz for mass, width, and coupling 
constant: 

M = Mz=\l + [Tz/Mzf] Mz^ Mz-l^^ Mz-3A MeV, (3.11) 



1 r3 



1 + iTz/Mzr\ 'Tz-Tz- -j^ -Tz-l MeV, (3.12) 

= . (3.13) 

" l + iTz/Mz ^ ' 

A naive use of a constant width in (|3.6|) - (|3.8| ) would lead to a wrong determination of what 
has been introduced there to be the Z mass. In fact, one can put forward a completely different 
point of view [Q-[jl4| (see also references cited therein and [^-[^). There is no physical 



reason to consider (^)-( P^ ) as the final result of a perturbative calculation. After so many 
formal manipulations, including renormalization, one could consider the transformations (|3.9|) - 
( p.l2|) as part of the renormalization procedure. In doing so, one is in complete agreement with 
the ideas of S-matrix theory: Unstable particles are described by simple poles of the S-matrix 



^ In the Standard Model, this is the case if two conditions are fulfilled: (i) there are no opening new Z decay 
channels (production thresholds) near s = M^; (ii) radiative corrections to Q are practically independent of s 
in a region where Q essentially influences the cross sections. 
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in the complex energy plane whose location is defined by the particle's mass (real part) and 
its width (imaginary part). Such an approach automatically anticipates the propagator ICz{s) 
with mass (|3.11|) and width (|3.12|) . 

The default mass and width definition for the S-matrix branch of the zpl Tj^j^ package is 
( p.9|) -( ^TT2D , while (|3.6|) - (|3.8| ) should be used for all other branches. The final choice of the 
definition of the Z mass and width is left to the user (see fiag GAMS). 

For the t channel of Bhabha scattering, the above discussion regarding the propagators is 
also of some relevance. Of course, here the width of the resonance is absent. Furthermore, 
one must replace s by t and s' by t' in the propagator functions. 

rp 

For the corresponding Standard Model calculations, zpl Tj^n relies on the DIZET pack- 
age. The following parameters are passed to DIZET: 



a,as,G^,Mz,MH,'mf, 



(3.14) 



which returns M^y; the total and partial Z widths; the weak form factors, etc. Thus, one 
arrives at improved Born cross sections, which are used as building blocks of the QED for- 
mulae discussed in the foregoing sections; zpl does not calculate bare Born cross sections 
since definitions of the latter tend to be ambiguous. The weak mixing angle also will not be 
considered as a quantity of physical relevance (although one could do so), but will only be 
used for book-keeping of intermediate results. 

In the remainder of this section and in the next, we discuss the various assumptions 
regarding the functions lA{m,n; s) that have been implemented in zpl Tj^^. 



3.1 0{a) Corrections to Ar 



In the on-mass-shell renormalization scheme |39| that is used in zpl Tj^^i the weak mixing 



angle is defined uniquely through the gauge-boson masses: 



sin^e 



w 



sin^ Ow 



Mr 



(3.15) 



(3.16) 



1 - Ar 

In subroutine SEARCH of DIZET, Ar is calculated to order 0{a) as defined in ||3^, where the 
heavy top contribution is calculated as in M. Recently, a careful comparison [HOl of two 



independent 0{a) calculations of (|3.16| ) showed agreement in 12 digits 



3.2 0(a) Corrections to Tz 



Electroweak corrections to the Z width have been calculated to order 0{a) in The partial 
decay width of the Z into fermions of type / is given by: 



G 

^Y^/i^QEDC//?QCD(M|)pf X 

I [l - 4\Qf\ sin^ ewf^f + S{\Qf\ sin^ Owt^ff 



V^247r 



/i -RqED Cf -Rqcd {Mz) 



+ a 



'Ml, 



f 



1 + 2 



f 

Ml, 



6(a 



Ml 



(3.17) 
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The renormahzed vector and axial- vector couplings^ are defined as follows: 



(3.1^ 



1 -4|Q/|sin2 Owti] 



(3.19) 



The bare Born vector and axial- vector couplings, v and a, correspond to p = k = 1. The 
weak form factors, p and k, are real, constant, and depend slightly on the decay channnel. 
They contain the electroweak corrections to the process, including the dependence on rrit and 
Mh- Since zfI^Tjt^^ exactly follows 0, we need not go into details here. However, it must be 
mentioned that the Z width depends on the choice of the definition of the Z mass (see the 
discussion presented at the beginning of this section). 
Sometimes, the combination 



2,efr 



sin^e 



w 



(3.20) 



is used as a definition of the 'effective' weak mixing angle; see e.g. [41, 42] and the recent 



discussion on possible alternatives in and references quoted therein. Such a definition can 
also rely on the other decay channels: 



2,f 



sin 6w 



(3.21) 



Several factors contain additional corrections: 



/i = /i(M|), 



1 



4mj/s, 



(3.22) 



R 



QED 



QCD 



1 + Ciiml 



3 a 
Ait 

[ 3 for quarks 
] 1 for leptons 



1 + 0.0017Q2 



+ 



TT 



for quarks 
for leptons. 



(3.23) 



(3.24) 



(3.25) 



The corrections -Rqed and -Rqcd contain the photonic and gluonic corrections, respectively. 
The factor ci is of relevance only for the b-quark channel. The exact definition of -Rqcd will 



be given in section 



Owing to the large mass splitting between the t and b quark, there are two vertex diagrams 
for the Z decay into b quarks (fig. 0), which contribute additional mj-dependent corrections 
that are absent in the cases of light quarks [0, ^ (see also 0). The matrix element 
for the decay of the Z boson into d and s quarks (and similarly, for u and c quarks) may be 
written as follows: 



Mr. 



u. 



(3.26) 



^The Born axial- vector coupling is often defined to be equal to the weak isospin, a/ = I^if), e-g- de 
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Figure 7: Top quark exchange diagrams which contribute to Tt,. 



sin^ 6w 






















Ttot 


0.2282 


166.6 


83.6 


83.6 


83.4 


296.6 


382.9 


296.2 


382.9 


0. 


375.7 


2484.7 



Table 1: Weak mixing angle and partial and total Z widths; widths are given in MeV. 



The corresponding matrix element for b quarks has an additional left-handed contribution: 

7/3^^f + 7/375^1 + A6("^?)7/3(l + 75) 



M, 



u 



~ ^ ^Mlt^^a.u [7/3(1 + 75) - 4 sin^ ^^/cf 7^3' 
Here, Aj, vanishes for m-f — > 0. By simple algebra, one can show that 

Pf - 2- 



u. 



(3.27) 



K 



1^^ + 



(3.28) 



(3.29) 



These exact form factors have been implemented in zpl Tj^^i- In the limit of large t-quark 
masses, the leading terms are given by 



ab 



a 



An sin^ 6 



w, 



tb\ 



w 



+ - + 



3 6 cos^ 6 



w 



log 



mi 



(3.30) 



rp 

where Vtb is the (t, h) Kobayashi-Maskawa mixing matrix element. zfI Tj^ji calculations use 
the normalization = 1 in ( p. 271) . 

The calculation of the W width follows the same principles as that of the Z width and 
is realized in subroutine ZWRATE of DIZET. Since the W width is not that important for the 
description of fermion pair production, we do not go into details. 
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The partial and total widths, returned by zfI Tj^jjl, are summarized in table |I|. In addition 
the weak mixing angle is given in the table. 



3.3 0{a) Corrections to Fermion Pair Production 

Fermion pair production in the Standard Model gets contributions from self-energy insertions, 
vertex corrections and box diagrams. We divide these into two gauge-invariant subsets. 

Fermion loop insertions to the photon propagator are summed together with the photonic 
Born diagram (see fig. ^) to form the matrix element Ai-y. In effect these corrections change 
a into a{s): 

1 



Mj{s) 



where the following short notation for bilinear combinations of spinors uj is used: 

Ai3®B'^= [UeApUe] " \ufB^Uf . 



(3.31) 



(3.32) 



After a Dyson summation of the fermion loop insertions Aa(s) to the photon self energy, the 
running electromagnetic coupling constant contains higher-order corrections: 



a[s] 



FA(s)a = 



a 



1 - Aa(s)' 



(3.33) 



The function XFOTFl in DIZET is used to calculate Aa. 

Some numerical examples are given in table |^ as a function of y/s. 





30 


87 


89 


91 


93 


95 


200 


Fa 


1.0504 


1.0630 


1.0633 


1.0635 


1.0638 


1.0640 


1.0723 




- i.0186 


- i.0188 


- i.0188 


- i.0188 


- i.0188 


- i.0189 


- i.0191 



Table 2: The running QED coupling, Fa{s) = a{s)/a, as a function of the centre- of-mass 
energy (y/s in GeV). 




7 





Figure 8: Photon (a) and Z (b) self energies. In (b), the case 7,7 is not included. 



^In the examples we have taken Mz — 91.175 GeV, Mh ~ 300 GeV, nit 
flag values (see table ^in section 3.2) unless explicitly stated otherwise. 



140 GeV, as = 0.12 and default 
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At ^/s = Mz the running coupling constant has the value: 



(3.34) 



In addition to the running of a{s), there are diagrams with additional internal W and Z boson 
propagators to the photonic Born amplitude, e.g. a ly^-pair insertion or a vertex correction 
with a Z propagator. These diagrams could be treated as corrections to the photon amplitude 
as well. However, this would make Ai^ dependent on the gauge. Diagrams of this type form a 
gauge-invariant subset together with all the insertions to the Z Born diagram as well as with 
ZZ and WW boxes. So, any diagram with at least one additional massive gauge boson will 
be combined with the Z exchange Born diagram to form the matrix element A4z (see figs. H), 



and |T0]). 




1,Z 





Figure 9: Vertex corrections to the Z matrix element. 



W 



W 



(b) 



Figure 10: Box diagrams contributing to the Z matrix element. 



The contributions to the corresponding matrix element Aiz can be expressed in terms of 
four weak form factors (pe/, i^e, i^f, t^ef) as introduced to order 0{a) in |P]: 



Mz{s, cost?) 



G f,aea fPefjs, cos ^) 
s- Ml + iMzTz 



Lp® L'^ - 4|Qe| sin^ 6wi^e{s, cos^)'^^ (g) L'^ 



-4:\Qf \ sin^ 6wKf{s, cos?9)L/3 7^ + lQ\QeQf \ sin^ 6wK,ef{s, cosi!})^i3 ® 7^ 



(3.35) 



Lf3 = 7/3(1 + 75), 



(3.36) 



where Lp is the left-handed projector. The (p — 1) and {n — 1) are normalized with the factor 
a /{in sin^ ^w)- 
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The matrix element may be rewritten in terms of renormalized vector (v) and axial- vector 
(a) couplings: 



Mz{s, cos-i?) 



G, 



s-Mi + iMzTz 



+aet;/7/375 ® 7^ + ® 7'^], 



(3.37) 



a/ = y^Pe/(s,cos^9) /3^(/), 



(3.38) 



1 — 4|Qj| sin^ 6wf^f{s, cos'i?) 



(3.39) 



fe/ = CLeVf + fefl/ — ttettf 1 — 16|(5e<5/| sin"^ 9wf^ef{s, COS ??) 



(3.40) 



The four form factors are the most general ansatz for the weak radiative corrections. In the 
Born approximation, p = n = 1, and Vef = VgVf. The coupling Vej has no parallel in the Born 
approximation and is, in principle, completely independent of Ve and Vf. Form factors are 
calculated in subroutine ROKANC of the DIZET package. 

In table ^ we show the s dependence of the weak form factors for lepton production. 





30 


87 


89 


91 


93 


95 


200 


Pef 


0.9992 


1.0020 


1.0021 


1.0022 


1.0022 


1.0023 


1.0102 


- i.0006 


- i.0043 


- i.0045 


- i.0047 


- i.0048 


- i.0050 


- i.0283 


Ke, Kf 


1.0283 


1.0227 


1.0226 


1.0226 


1.0225 


1.0225 


1.0117 


+ i.0115 


+ i.0134 


+ i.0135 


+ i.0135 


+ i.0136 


+ i.0137 


+ i.0329 




1.0552 


1.0459 


1.0458 


1.0456 


1.0455 


1.0454 


1.0355 




+ i.0204 


+ i.0265 


+ i.0268 


+ i.0271 


+ i.0273 


+ i.0276 


+ i.0560 


KeKf 


0.0020 


-0.0002 


-0.0002 


-0.0002 


-0.0002 


-0.0001 


-0.0010 




+ i.0032 


+ i.0008 


+ i.0007 


+ i.0006 


+ i.0005 


+ i.0004 


+ i.OllO 



Table 3: Leptonic form factors pef, He 
m GeV). 



tif,Hef as functions of the centre- of-mass energy (y/s 



The last row in the table shows how well Kg/ can be factorized in terms of Kg and Hf. 



The 

form factors shown in the table have been calculated without the box diagrams of fig. [1^. The 
program DIZET allows for three options: inclusion of these box diagrams into the weak form 
factors; calculation of them as an extra cross-section piece, o"box, to be added incoherently (see 
[^); or neglecting them completely. zfI Tj^j^ allows for the last two options only. This has 
been arranged in order to make the weak form factors independent of the scattering angle. 
Thus, the angular integration could be performed analytically. Table § shows the influence 
of the box-diagram corrections on the differential cross section for different angular bins at 
LEP I energies^. 

Henceforth, we omit the possibility of an angular dependence of the form factors. Using 
an alternative parametrization, the axial-vector couplings may be chosen such that they are 

^ The contributions from box diagrams are non-resonant at LEP I energies. zpf^Tj^j^ users should be aware 
that, off the Z resonance peak, box diagrams may not be neglected with respect to other radiative corrections. 
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1? region 




= 94 GeV 


O'fM — Chox 




rr, v 1 n''' 
(^box X -LU 


O'n — CThox 






u — ou 


0.1410419 


0.1410410 


9 


0.0512613 


0.0512605 


8 


30° - 60° 


0.2993644 


0.2993631 


13 


0.1042011 


0.1041998 


13 


60° - 90° 


0.3005794 


0.3005789 


5 


0.0961092 


0.0961088 


4 


90° - 120° 


0.3016901 


0.3016900 


1 


0.0800074 


0.0800073 


1 


120° - 150° 


0.3015250 


0.3015249 


1 


0.0721641 


0.0721640 


1 


150° - 180° 


0.1420456 


0.1420456 





0.0353764 


0.0353763 


1 



Table 4: Differential cross section in nb for fi^fi production both with and without weak box 



contributions for selected LEP I energies. QED corrections are calculated following section 
with = 4mj/s. 



unchanged by radiative corrections. In this case, the Fermi constant absorbs the weak form 
factor p{s) and becomes dependent on the process and its kinematics: 



— ^ — p{s)G^. 

The other form factors may be absorbed into various weak mixing angles: 

Ke{s) sin^ 9w 
sin" ew = l- M'^/Ml { ^ /(g)sin ^ 

'Kef{s) sin^ 6w 



(3.41) 



(3.42) 



This is similar to ( ^.2U| ), even though more involved because of the additional complications 
presented by the kinematics. The above parametrization of Aiz allows for a Born-like in- 
terpretation of all weak corrections. In this respect, we differ in our intentions from many 
other definitions of weak form factors and couplings, which try to perform dedicated approx- 
imations. Of course, such approximations may be applied to our weak form factors or to 
quantities derived from them; in sections ^ and |8]^ such approximations will be discussed; see 



also section p.2| . At LEP I energies, the approximate relations hold: 

Pe/(M|; 

where the second relation may be replaced by: 



~ K 



f 



(3.43) 
(3.44) 

(3.45) 



So far, we have concentrated on s-channel kinematics, which depend on s and cos-i?. It 
should be noted that for the t channel in Bhabha scattering the energy variable becomes 
t = — |s(l — cost?) instead of s. 

For b-quark production, unlike d- and s-quark production, a special contribution to the 
weak form factors arises from diagrams in fig. Pd, which contain as building blocks the Feynman 
graphs of fig. |^. This contribution may be of special interest at a high luminosity version of 
LEP I as is discussed in [|^. In general, the correction is s-dependent. It can be approximated 
near the Z resonance by the corresponding correction Ab{mf) to the Z width as introduced 
in (p^): 



A 



Peb — Ped 



bjjnt 
ah 



(3.46) 
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Kb = Kd+ 



ab 
ab 



(3.47) 



(3.48) 



with Ke unchanged. This approximation has been implemented in zfI Tj^j^', it is vahd only 
near the Z resonance and for rrit > y/s/2. At other energies, since it would be difficult or 
impossible to measure the effects of this tiny correction due to small cross sections, we assume 
that the approximation holds there as well. 

Higher order corrections in -Rqcd have been implemented in zpl Tj^j^- These and the effects 
of the higher-order corrections to p and k associated with a potentially large t-quark mass will 
be discussed in the next section. 

We now come to the cross-section formulae, which are calculated in subroutine BORN. Both 
(Tt and apB are sums of three terms: 



o,SM/ \ 



0,SM/ \ , 0,SM/ ry\ 



iri7,7;s)+^e irh,Z;s)}CUs) + iTiZ, Z; s) \}Czis)\\ (3.49) 



rSM 



0,SM/ ry r7\ 
a A [S] Z, Z) 

rSM/ 



Here, we have written the cross section in a form that is suitable for initial- and final-state 
radiation. However, the initial-final state interference cross section depends on two different 
energy scales (s, s'). The correct handling of the propagators can be inferred from ( |3.1| ). The 
generalized couphngs X4 are assumed to be dependent on s in zpl Tj^^ with the exclusion of 
the running QED coupling where the scale (s or s') can be chosen by a flagj^. This assumption 
speeds up the calculations with negligible loss of accuracy. In principle, one can take into 
account the s and s' dependence, in a trivial way for the factorizing parts of the form factors, 
and the rest with a little effort. 

For unpolarized scattering, ctt can be expressed by 



(3.50) 



J|^(7, Z- s) = 2cmN^{s)Nz\QeQf\ [FXp 



e} Ve}\ , 



(3.51) 



Xf^{Z,Z-s) = N,{s)Nl 
where 



/i(s)c/ i?QCD(s), 



v/2 87ra' 
;i + 2mj/s). 



\Pef\ 



(3.52) 

(3.53) 

(3.54) 
(3.55) 



The variables /i(s), c/, and -Rqcd have already been introduced in section |3^ , and F4 in 



^ In principle, with initial-state radiation, the form factors depend on s', with final-state radiation on s and 
with the small initial-final interference on both s and s' . 
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The corresponding generalized couplings for the anti-symmetric cross section apB are: 

lfi{l,r,s) = 0, (3.56) 
lf^{j,Z;s) = 2fi{s)N,{s)Nz\QeQf\PefFX, (3.57) 

ifiiZ, Z- S) = Ms)N^{s)Nl{VeV} + Vef)\pef\'. (3.58) 

In apB, the QCD-factor is set zero, -Rqcd = 0. 

Helicities and polarizations may be included in the Standard Model cross sections cta in a 
compact way for massless fermion production P, |[. To do this, one must replace the above 
couplings, TA{m, n; s), with|^ 

Cr(m,n; Ai, A2,/ii,/i2) = 

{Xi[ve{m)v*{n) + ae{m)al{n)] + X2[ve{m)al{n) + vl{n)ae{m)]} x 
{hi[vf{m)vj{n) + af{m)a}{n)] + h2[vf{m)a*f{n) + Vj{n)af{m)]}, (3.59) 

CpBim, n; Ai, A2, = Crim, n; A2, Ai, h2, hi). (3.60) 

The vector and axial- vector couplings Vf{0) and a/(0) of the fermion to the photon are: 

VfiO) = QfFA{s), a/(0) = 0. (3.61) 

Here, we introduce the longitudinal polarizations of the electron (A_) and positron (A+) 
and the helicities of the final state fermions h± in the following combinations: 

Ai = l-A+A_, A2 = A+-A_, (3.62) 

hi = ^{l-h+h^), h2 = ^{h+-h_). (3.63) 
The various parts of the cross section in { ^A\ , p.49|) now become: 

cT^(s, s'; m, n) = 3?e |cA(m, n; Ai, A2, h, /i2)^[Xm(s')Xn(s) + Xm(s)x* (s')]| , (3.64) 



X-y{s) = yiV^/C^(s), (3.65) 

Xz{s) = ^N^{s)NzlCz{s), (3.66) 

where again the QCD-factor -Rqcd in is set equal zero for apB- 

If at least one incoming and one outgoing fermion are polarized, then the contribution to the 
forward-backward anti-symmetric Standard Model cross section from pure photon exchange 
does not vanish as in ( p.56| ). This can be seen from formulae ( p.59| )- p.64| ). 

We have not gone into details that are specific to Bhabha scattering: this is done in [|ri|] 
and it represents a straight-forward extension of the above formulae. 

This is not rigorous with respect to Vef, which has been assumed to factorize in order to simpHfy the 
notation. The correct expression, implemented in zpl Tj^j^, can be obtained by performing the multipHcations 
in ( 3.59| ) and re placi ng the product Ve{Z)vf{Z) with Vf,f{Z). This may be verified by explicitly squaring the 
matrix element ( |3.35 ). 
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Asymmetries represent a clean and near systematic free measurement with which to test 
various models. In addition to the forward-backward asymmetry, several other asymmetries 
are interesting. It is useful to define a generic 'spin' asymmetry, A{h): 



A{h) 



a{h)-a{-h) 



a{h) 



-hV 



(3.67) 



where h can denote the polarization of an incoming fermion or the helicity of an outgoing one. 

Choosing h to be /i+ = +1 the helicity of a final-state r"*" and to be aT{h+), A{h) 
becomes the r polarization, = Apoi- Similarly, one can define: A^^ as in |^ from 

( p.67|) with a{h) = aA{h+),A = F,B,FB, respectively. The subscript F (B) is used to 
indicate that only data from the forward (backward) hemisphere are in the measurement; 
the corresponding theoretical relations are given by ( |3.64| ) and (|1.6| ). The forward-backward 
polarization asymmetry may be defined as follows: 



^FB — ^ 



(3.68) 



3.4 Higher-Order Corrections 

Here, we give a summary of treatment and common resummation of some higher-order weak 
and QCD corrections in zfI Tj^^. 

Some higher-order terms are used to correct Ar, p and k. These terms are exclusively due to 
t-quark mass corrections. zfI Tj^j^ takes into account the following mj-dependent corrections: 

• complete m^-dependent 0{a) terms 0, 

• leading 0{a^m^) terms [||, 



complete (either approximated as a function of energy or exact) 0{aas) terms |]SD|, ^ 
with leading part O^aagm^) . 



For Ar as introduced in (|3.16| ), a common resummation of these leading terms may be 
performed as follows |5|, |0|, ||0: 



1 



Ar 



cos^ 9 



[1 - Aa(M|)] (1 + ^-^^6p) - Ar, 
sm uw 



Ar- 



^^lloop 



cos^ 6 



w 



rem -2/1 

sm (J 



w 



Ap - Aa(M|) + Ar; 



etas 
rem ; 



(3.69) 



(3.70) 



Ap 



Ap° + Ap"^ 
3a 



Xn 



rat 



167r sin^ Ow cos^ Ow Ml 



1 _ -(1 I ^% "^(g^'^Ms) 



2 

3' 



TT 



(3.71) 



A detailed discussion of the zpI^Tj^ji flags which control the implementation of these corrections will be 
presented in section 6.2. 
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where the Aa{s) is introduced in ( ^.331 ) and 

Xo = 



Ml 



Hoop 



(3.72) 



MS 



In Xq the UV divergencies are removed according to the MS renormahzation scheme with 
/i = Mz Q The separation of Xq is not uniquely defined; it introduces a dependence of the 
resummation on the renormahzation procedure. Further, 



6p 



3r 



1 - {2n^ - 19)r- -(1 + 



r 



X = 2 sin^ 9w cos^ 9w 



3' 3 
v^87r2' 



TT 



1 - Aa(M|) 



(3.73) 
(3.74) 
(3.75) 



For the cross-section form factors, p(s, cos??) and k(s, cost?), and partial Z width form factors, 
p^ and K'^, similar formulae hold: 



P 



„lloop I aois 
r ' r rem 



Ap 



lloop _|_ aas 



1-Sp 

cos^ 9w 
sin^ e 



w 



Ap] 1 



cos^ 6w 
sm 6w I 



For the cross section alone, we have additionally: 



W 



Ap 



^1 + '^s-p 



sin^ 61 



(3.76) 
(3.77) 

(3.78) 



Some numerical examples of the effect of leading 0{a'^mf) corrections on the weak mixing 
angle, the p^ p~ total cross section and forward-backward asymmetry at s = M| are shown 
in table |^. 

In section we introduced the QCD correction factor, -Rqcd, in ( |3.25| ). Its exact definition 



as implemented in zpl Tj^^ is given by [ p5| : 



QCD 



1 + — + 1.409 



TT 



12.805 



a. 



TT 



QCD3. 



(3.79) 



For b quarks, the top- and bottom-quark mass- dependent QCD corrections (ci,C2) up to 
0{a'^) have been taken from 



a. 



-Rqcd = 1 + ci(mfe) h C2(mb, rrit) — 



TT 



a. 



TT 



12.805(^1 QCD3. 



(3.80) 



Where QCD3 has the value or 1 as required by the user. In subroutine ZUWEAK, -Rqcd is 
calculated and the result is stored in the variables QCDCOR and QCDCOB. 
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In |B9[ this corresponds to a replacement of Mw by Mz in the UV divergence T'uv- 
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\J UocI VcLUlc 




(GeV) 
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^ /1 8/1 7 
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9 


1.4847 


1.4867 


1.4895 


1.4935 


1.4990 




3 


1.4847 


1.4867 


1.4896 


1.4935 


1.4990 







-.0043 


-.0024 


-.0004 


.0015 


.0023 


^FB 


1 


-.0045 


-.0024 


.0004 


.0040 


.0078 




2 


-.0045 


-.0023 


.0008 


.0052 


.0115 




3 


-.0045 


-.0024 


.0007 


.0050 


.0113 



Table 5: The weak mixing angle, muon pair production cross section and asymmetry both with 
(AWT^^O) and without fAMT4=Cj leading 0{a^mj) terms. 



Depending on a flag, the running strong interaction coupling constant |j55| a(g^,Aj^) is 
calculated with functions ALPHA4 or ALPHAS [El : 



An 



1 - 



\nA + 




(3.81) 



with 



bo = 11- -nf, 



6i = 102-yn;, 



„ 5033 

2857 Uf 

9 ^ 



325 2 
27 ^. 



A = \n- 



{n;)2 ■ 



MS 



(3.82) 
(3.83) 



Here represents the energy scale and rtj the number of quark flavors. The corresponding 
definition of may be found in table |[ 

In -Rqcd, the is calculated with = s for cross sections or = M| for the Z width. 
In the 0{aas) corrections the scale is chosen to be = M| for light quarks and = 
max{M|,mj} for the tb doublet. One can see that in the LEP energy region the difference 
between the approximate]^ and exact treatment of the 0{aas) corrections is minor. 

Table shows the dsil^, ^Ms) ^ function of g^. 

The dependence of weak parameters on is shown in table |^ as a function of how the 
0{aas) corrections are apphed. 
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The approximation realized in ZF^Tj^ji is a Taylor series expansion in s/rri^ for s smaller than m^. The 



leading term of this expansion is included in ( |3.73 ) 
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< mi 



< mi 



> mi 



MS 



A (4) ( 
MS I Zll 

\ MS J 



2_ 

27 



107 
2025 



MS \ nif, 



In 



MS 



(4) \ 23 
MS 1 In 



> MS' y 



MS I mt 



(5^\ 2T 
MS 1 \^ 



(a1£L)2 



-107 
1127 



Table 6: Ky^^ for different energy regions. 



\q\ (GeV) 


50 


100 


150 


200 


250 


300 




.1193 


.1079 


.1022 


.0985 


.0958 


.0937 



Table 7: Running as{(f ) versus \q\. 
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QCDC 










vv final state 





1.004 


- 10.002 


1.025 + 10.013 






1 


1.004 


- 10.002 


1.023 + 10.014 






2 


1.004 


- 10.002 


1.022 + 10.014 






11 final state 





1.003 


- 10.005 


1.025 + 10.013 


1.025 + 10.013 


1.050 + i0.026 


1 


1.002 


- 10.005 


1.023 + 10.014 


1.023 + 10.014 


1.046 + 10.027 


2 


1.002 


- 10.005 


1.022 + 10.014 


1.022 + 10.014 


1.045 + 10.027 


uu final state 





1.003 


- 10.004 


1.025 + 10.013 


1.024 + 10.012 


1.050 + 10.026 


1 


1.003 


- 10.004 


1.023 + 10.014 


1.022 + 10.013 


1.045 + 10.026 


2 


1.003 


- 10.004 


1.022 + 10.014 


1.022 + 10.013 


1.045 + i0.026 


dd final state 





1.003 


- 10.003 


1.025 + 10.013 


1.024 + 10.012 


1.049 + 10.025 


1 


1.003 


- 10.003 


1.023 + 10.014 


1.022 + 10.012 


1.045 + 10.026 


2 


1.003 


- i0.003 


1.022 + 10.014 


1.022 + 10.012 


1.044 + 10.026 


hh final state 





0.999 


- 10.003 


1.025 + 10.013 


1.028 + 10.012 


1.054 + 10.025 


1 


0.999 


- 10.003 


1.023 + 10.014 


1.026 + 10.012 


1.049 + 10.026 


2 


0.999 


- 10.003 


1.022 + 10.014 


1.026 + 10.012 


1.049 + 10.026 



Table 8: The dependence of the form factors Pe/ o.nd k^, kj, K^f at the Z peak on different treat- 
ments of 0{aas) corrections (13CDC=0 - no, 1 - approximate and 2 - exact 0{aas) corrections 
are applied). 



30 



4 The Hard-Scattering Process: 
(II) Model-Independent Branches 

4.1 Effective Couplings 

In a simple quantum mechanical approach, the Z boson may be assumed to have real constant 
vector [vf) and axial- vector (a/) couplings to fermions (/). This ansatz may be realized by a 
replacement of the renormalized effective couplings as predicted from the Standard Model by 
naive effective couplings in the cross section, see fig. |TT]. 




Z 




Figure 11: Matrix element in the approach with effective Z couplings. 



The cross sections for the hard-scattering process are: 



o.cS I \ 



0", 



\s] 7, 7) + CT°f^{s] 7, Z) + a°/^{s] Z, Z) 



= (7, 7; s) + 3?e [Xf (7, Z- s))CUs)] + ifiZ, Z; s) \ICz{s) 
Here, the generalized couplings for the total cross section are: 



lf{j,r,s) = c^N,{s)\FA{s)fQlQ' 



(7,Z;S) = 2c^NzN^{s)Fa{s) \QeQf\VeVf, 



lf{Z,Z-s) = NlN,{s) 



Cm{vl + al){v) + a)) 



where N^{s) and Nz are defined in ( |3.53| )-( |375^ ), Cm in ( p.55| ), and Fa in ( p.33| ). 
The asymmetric cross-section part is defined by: 



If^i'y, Z; s) = 2fi{s)NzN^is)FAis) \Q,Qf\ a.aj, 
Tf^{Z, Z] s) = 8n{s)NzN^{s) VeaeVfaj, 



(4.1) 

(4.2) 
(4.3) 
(4.4) 



(4.5) 

(4.6) 
(4.7) 



where /x(s) is defined in ( |3.22| ). 

One may interpret effective couplings as approximations to the weak form factors of the 
Standard Model e.g.: 

' (4.8) 
1)1 , (4.9) 



3?e 
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where the alternate notation {(][,(]:[) is favored by the LEP experiments. We neglect here 
possible dependences on the scattering angle. In zfI Tj^ji the normalization a/=l is used for 
all fermions. In addition, one may choose an alternative parametrization in terms of the 
effective weak neutral current amplitude normalization, pj: 



Pf 



9 a 



PePf 



Pef 



(4.10) 



(4.11) 
(4.12) 



In the present approach, one can leave T z as a free fundamental parameter of the ansatz. 
Alternatively, one may define it through the second line of ( p.l7|) , replacing there the renor- 



malized (f , a ) couplings by effective 



') ones. 



In either case, one must realize that the normalization of the Z width may change de- 
pending on the definition of Mz (see earlier discussion in connection with equations (|3.8|) 
and ( |3.12| )). For additional general comments on cross sections and asymmetries we refer to 
section |^. 

In principle, this branch is completely model independent. However, zpPTj^fi users should 
be aware that the current implementation contains small Standard Model contributions in the 
form of imaginary parts of weak form factors (see table |^). 



4.2 Partial Z Widths 

This approach to determining cross sections relies on the assumption that scattering through 



the Z boson may be considered as subsequent formation and decay of a resonance, see fig. 12 



The corresponding net cross section is (see e.g. and references therein): 

a"/{s) = a^^^{s; 7, 7) + /"^ + cr^'^(s; Z, Z), 
(T^^{s; Z, Z) = lJp{Z, Z] s) |/Cz(s)p, 

The photonic contribution to the cross section, cr^'^^, is given in ( p.50|) - (|3.51 



(4.13) 
(4.14) 
(4.15) 





Figure 12: Matrix element for the Z resonance scattering. 
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One of the attractive features of the partial widths approach is the treatment of hadron 
production. Here one simply replaces Tf with Fhad in ( [4.15|) . In case of hadron production, 
we have to use the quark language for the photonic cross-section contribution in ( [4.13|) , which 
is taken over from ( |3.49| ) in section without changes. 

The interference term in (^4.13|) presents some complications with this approach. There are 
at least four different ways dealing with this term: an exact calculation, a simple parametriza- 
tion, ignoring the term completely, or assuming the Standard Model prediction. 

A correct handling would rely on partial decay widths into specific helicity states, r-|-(/). 



as proposed in |T2 



X^(7,Z;s) 



= (^T (s; 7, z) = 

3 A7rQeQfa{s)c„ 



nif)-Tlif) 



(4.16) 
(4.17) 



"8s MzRqed 

In practice, however, experimental measurements seem not to deliver a sufficiently high number 
of degrees of freedom to use this formula. 

A simple parametrization of the interference term, as has been realized in the S-matrix 
approach, could be used. This, however, leads to large uncertainties in the Z mass determi- 
nation [jl4| . 

Ignoring the interference term completely ( i.e. assuming that it is identically zero) is 
another possibility, since this term is expected to be small. In addition this removes one 
degree of freedom. 

The last alternative is to assume the interference term from the Standard Model, as defined 
in (H), (pTJ): 



o,SM, 



(4.18) 



This is the approach that has been implemented in zpl Tjr;^. 

In addition to the partial decay widths of the Z, the total Z width and mass are free 
parameters with this approach. There is an ambiguity in the normalization of the total and 
partial widths which is due to the different choices for the definition of the Z propagator 
(see discussion at the beginning of section The energy-dependent total Z width, r^(s), is 
related to the constant total Z width, Tz, by (p.8|), ( p.l2| ). In relating the two approaches to 
the resonance definition, there is no explicit constraint on the partial widths. 

They may be related as follows. In the Standard Model branch, the residua of the resonance 
functions were normalized by the Fermi constant. With this approach the actual residua are 
contained in the partial widths. Comparing the two, one can derive the relation between the 
partial widths for the two different definitions of the Z propagator from (|3.13|) : 



1 + {Tz/MzY 



(4.19) 



This relation is in full accordance with the corresponding relation for the total width ( 3.12|) 



4.3 S-Matrix 

Besides the approach to the effective Born cross section based on the Standard Model or on one 
of its extensions, there is only one accurate model-independent approach to the Z line shape. 
One can derive this rigorous model-independent formula either starting from an analysis of 
the Standard Model results or from S-matrix theory. 
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An early application of the S-matrix formalism to LEP I physics may be found in . 



Recently, it has been proposed to use S-matrix theory for a global description of the hard- 
scattering process |3^, Such an ansatz has the advantage that it contains no special 
assumptions on the dynamics beyond general principles and the existence of both photon 
and Z boson. In the necessary formalism has been described. One starts from the 



incoherent sum of four squared matrix elements for the scattering of helicity fermion states 
{^L^\ JlIr, ele^ fafi, e^el /ij/i, e^et fLh) as seen in fig. 




Figure 13: Scattering in the S-matrix approach. 
All have the following structure: 

Mi{s) = ^ + + F.isz) + {s~ sz)Fl{sz) + {s- szfFl\sz) + • • • (4.20) 
s s - Sz 

The location of the Z-boson pole is given by sz in the complex energy plane, and and i?^ 
represent constant complex residuals. The cross section is: 

^f{s) = -j:s\M\s)\\ (4.21) 
^ 1=1 

In order to fit the cross section, it is perhaps useful to decompose the above into a series of 
real-valued terms with rising powers of [s — M|) |^: 

, M%R + is — M'l) I rn , os. r-i , o^o To 

= f + |.-,,|^ + 4 ^ ^ "'^jk ^ ^ ^'^^M ^ ■ ■ ■ 

where we defined the constants to be dimensionless and 

= (4.23) 



f> ^, Qe \/3Qd + 2gVf c/^QCD «(M|) for hadrons at LEP I 
^ QeQf^a{Ml) forleptons. 

Here depends only on the dynamics of the photon, while the other parameters (M^, F^, i?, /, 
ro, . . .) depend also on the Z. The ansatz ( [4.22|) may be compared with the notation used in 
the foregoing sections: 

a°/{s) = a^'^(s;7,7) + a^'^(s;7, Z) + (Ty'^(s;Z, Z) 

= J|(7, 7; s) + [j|(7, Z; s)}Czis)] + I^Z, Z; s) \ICzis)f . (4.25) 
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Each of the coefficients itself is then a series in the variable Tz/Mz- 
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Here, 



j|(7,Z;s) 



s 



J I-R 



+ 



s 
R 



+ 



(4.26) 
(4.27) 
(4.28) 



The dots indicate contributions from the Taylor coefficients; for instance, the photon-exchange 
contribution collects not only the but also small additional terms due to the dependence 
of the running QED coupling on s or s' The leading contribution to the jZ interference 
comes from the combination J = I — R\ m order to simplify the s dependence of the ansatz, 
R has been introduced in ( |4.22| ) with a coefficient of M| instead of s. 

In order to get an intuitive feeling for the meaning of the S-matrix parameters, it may be 



helpful to contrast this approach with that of the effective couplings discussed in section |4?1. 
The cross section for muon production (|4.22| ) in this approximation is given by fixing: 



R = c^{vl + al){vl + al), 



I = R + J, 



J 



2cR^VeV^,, 



0, i 



0,1,... 



with 



/47r G^Ml 



3 V2 Stt ■ 

Contrasting the S-matrix to the partial width approach, one obtains instead: 

R = 127rrer/ + . . . , 



(4.29) 

(4.30) 
(4.31) 
(4.32) 

(4.33) 
(4.34) 



and the 7Z interference part has to be fixed by a relation analogous to ( [4.30D (see the lengthy 



discussion in section [4.2| ). An exact treatment would enlarge the number of parameters to be 
fitted, and weaken the numerical result. 

The general form of the above parameters [R, /, tq, . . .) may be found in |1^ (they were not 
made dimensionless there as is done here). In a quantum field theory, the constants (ro, ri, . . .) 
are non-vanishing, owing due to non-resonating quantum corrections. A careful analysis of 
their calculation in accordance with the S-matrix properties has been performed in [p5[] . 

An ansatz quite similar to ( |4.22| ) has been derived in starting from an on- mass-shell 
renormalization of the Standard Model; for the production of flavor /: 



0-7 



o,S, 



\s - Sz[ 



Ml 



+ 7^ 



s-Ml 



f- 



+ 



(J, 



QED 5 



(4.35) 



where terms of higher order in (s — M|)/M| and in Tz/Mz are dropped. There is a simple 
one-to-one correspondence to the terms in ( [4. 221 ), with exclusion of the Xj; the dominating 



Strictly speaking, the residuum of the photon pole (4.24) is not a{M'^) but the QED coupling constant a 
at zero momentum; the difference is related to non- leading terms and may be absorbed by a redefinition of r^. 
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part of this correction is due to the imaginary part of the running QED couphng a{s). The 
corresponding contribution in our notations may be found in the exact definition of R: 



R 



W\R^z\' + 2^' 



(4.36) 



As in the aforementioned branches, the definitions of mass and width of the Z boson are 
correlated and deserve special attention. 



The possibility to describe asymmetries is mentioned in |jT4|. 



5 Beyond the Standard Model 

In recent years, many searches for possible effects from New Physics have been undertaken in 
precision experiments at LEP I. Rewievs of the present status and of the literature may be 
found in ^l[|-|6^. Here, we would like to restrict ourselves to some hints on the possible use 
of ZF^ Tj^ji for corresponding searches. 

To start with, let us assume that some more general theory leads to predictions for the 
scattering of two fermions into two fermions. This may be described by an additional matrix 
element J^e, to be added to and J^z- 



Me{s, C0S^9) ~ CE[UpLp ® L'^ + Me7/3 (^Lf^ + UfLp 7^^ + Me/7/3 ® 7^], 



(5.37) 



and the CE,Ua can depend on s and cos??. An instructive example for new physics at the 
Born level is an additional heavy neutral gauge boson Z' with mass Mz', width Tz', vector 
and axial- vector couphngs v'^.a'f, and coupling constant gz'. As long as this Z' does not mix 
with the ordinary Z, the influence on the scattering process is due io M.e and thus limited to 
higher energy 



Czi 



9% 



s - Ml,+iMz'Tz'' 



(5.38) 



Up = a'^a'fp'^f, Uf = {v'j - a'j)p'^fK'f, u^f = {v'^ - a'J{v'f - a'j)p^jK'^f, (5.39) 

with p' = = 1 for Z' Born physics. 

In general, the corrections Ua may also be due to some loop insertions to the Z matrix 
element from a generalized renormalizable theory, or even simply due to some Standard Model 
corrections not yet included in the definitions of the weak form factors. An example of the 
latter case had been given in the Z- vertex corrections from t-quark exchange, see (|3.46|) - (|3.48|) . 
Usually, the additional loop corrections are incorporated into the Z matrix element ( p.37|) : 



MiZ,E) = Mz + ME, 



(5.40) 



M{Z,E){s, cos^) ~ Cz 
+ 



(^eaf + — 



7/375 ® 7^75 + 



- - , I 



— i I ^ 



7/375 ® 7 + 



7/3 ® 7 75 
7;3® 74,(5.41) 
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(5.42) 



In case of loop corrections to the Z propagator, the ratio Ce/Cz is free of the resonating s 
dependence around the Z peak, and in some scenarios the corrections are even constant. 

In general, however, this is not the case. Coming back to the example of an additional 
Z\ there is evidently a potentially remarkable s dependence of the insertions, being even 
resonating near the Z' peak. 

In any case, one can go a step further and include the non-standard corrections into the 
form factors introduced in ( |3.35D by the following replacements: 



Pef Pef{Z,E) = pefil 



Ur 



Cz PefO'eO'f , 



pef 



Kef 



'"ef 



Pef{Z,E) 
Pef 



1 + 



PefiZ,E) 



^ef 



1 + 



Cz PefK.fae{Vf - 
Ce Uef 



«/), 



Cz pefl^efiVe - ae){Vf - ttf) 



(5.43) 



(5.44) 



(5.45) 



The above replacements ensure the interpretation of the weak form factors as finite renor- 
malizations of Fermi constant and weak mixing angle; see ( |3.41| ), ( ^.42] ). They can, however, 
drastically change the numerical behaviour of the form factors, which now need no longer 
be small. The advantage of the above formulae is two-fold. Besides the compact notation 
and simple interpretation, they may be used not only for the description of the fermion pair 
production process ( |1.1D . Without changes, they describe also the effects of new physics in 
Bhabha scattering or in the crossed channel, i.e. ep scattering. 

Besides the neutral current amplitude A4e, new physics may show up also in other phe- 
nomena, thus influencing fermion pair production in an indirect way. It is well-known that 
e.g. a Z', which mixes with the ordinary Z boson, may influence the Z and W mass ratio 
and the Z vector and axial-vector couplings - it is these effects, which may be searched for 
at LEP I. How they can be covered in the language of form factors has been explained in the 
references quoted above. In addition, a careful derivation of the weak form factors following 



the notations used in the present paper may be found in |T^, where the use of zf' for a 
Z' search is explained^. The main consequences are contained in the following replacements 
in the definitions of weak form factors , for partial widths and pef, k/, Hef for the cross 
sections: 



Pef 

l^f 

Kef 



Pmix(l - l//)Vf , 
(1 -S/)Kf, 

Pmix(l - ye)(l - yf)pef, 
{l-Xf)Kf, 

{l-Xe){l-Xf)Kef. 



(5.46) 



(5.47) 



Here, Xf,yf are small corrections to the Z-boson couplings due to the Z, Z' mixing, and pn 
is due to the related slight Z mass shift: 



Pmix 



M| cos2 9 



w 



Ml 
Ml 



M 



'Ml 
Ml 



(5.48) 
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See also the package ZEFIT M| at ZFITTERQCERNVM. 
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The parameter Mq is the Z mass of the standard model without Z' mixing. The pmix influ- 
ences the widths and cross sections directly, since we have replaced in ( p.l7| ) and in (|3.35|) , (|337|) 
the coupling constant a of the on- mass-shell scheme by the Fermi constant G^. These are 
related as follows (see ( |3.16| )): 



7ra 



2 sin Qy\/ cos^ Q 



w 



^Mo^(l-Ar) = ^M|p„ 



'1- Ar) 



(5.49) 



In the same way as (1 — Ar) becomes part of and pe/ without mixing, the factor pmix(l — Ar) 
becomes part of the form factors when Z and Z' mix. 

Another, completely different source of deviations from the standard model are self-energy 
corrections 11 due to new physics, which may lead to the following changes of the weak form 



factors [42 



Ape/(s)=Ap(0) + n^z(Mi)- 
Ap^(s) = HwwKM^) 



sTlzz{s) - MlTlzz{Ml) 



s-Ml 



Ak{s) = Ak(M: 



cos 6 



w 



sin 9 



w 



n^z(s) -n^z(M: 



(5.50) 
(5.51) 
(5.52) 



For completeness, the corrections for the charged-current form factor have been added. The 
corrections at LEP I may be obtained from the above expressions by setting s = M|. 

A different starting point has been used e.g. in ||6y, 0. There it is studied how one can 
disentangle new physics from the possibly large, unknown t-quark corrections of leading order 
Gfj^mf; see section |3.4 For this purpose one can introduce three new parameters: 



€2 = CgAp + 



ei = Ap, 



(4 - si) 
es = c^Ap + {cl - sl)Ak' 



(5.53) 

(5.54) 
(5.55) 



rp 

The quantities Ap, Ar^y, Ak', Sq may be identified with quantities used in zfI Tj^ji and intro- 
duced aboveFl: 

Arw = 1 - (1 - Ar)- ^ 



a 



9a 



VP 
2 



1 



'1 



1 



Ap), 



^ _ 1 4„2,cfr 

— -L 4:^^ 

ga 



,2 2 



2' 2 

1 -4(1 + A/t>, 



SnC, 



(5.56) 

(5.57) 
(5.58) 

(5.59) 



It is up to the user of zfI Tj^^ to decide which of the various coupling definitions available 
in the program and described in sections ^]2[ ^l3| , |^ are used as couplings Qv^Qa in the above 
definitions. The running QED coupling is defined in ( |3.33| ). Thus, while it may appear that 

''Here, we exactly follow the notations of i^]; see also the package ZFEPSLON ^ at ZFITTER@CERNVM. 
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the e parameters are merely rearrangements of previously defined quantities, their merit lies 
in separating out the mt-dependent effects in ei and other (Higgs) effects in €3. Furthermore, 
for an analysis of LEP I data alone, the €2 parameter may be ignored. 

Besides the notations introduced so far, there are several similar ones used in the litera- 
ture, often in quite a different context. As one important example, we quote the following 
notation which introduces again some self-energy corrections, but now calculated before 7, Z 
mixing ||63|| : 



aT 



sin^ OwM^ 



[nii(o)-n 



33 



all f» 4e 



d 



nii(g' 



n 



331 



The relation to the above notations is ||6^ 
ei = aT, €2 = 



" U 



a 



£3 



4 sin^ 9w 



-S. 



(5.60) 
(5.61) 
(5.62) 

(5.63) 



Further relations between different notations may be found in |61|. 



6 Initialization of zfI 

zpI^Tj^ji is coded in FORTRAN 77 and it has been implemented on IBM, IBM PC, VAX, and 
APOLLO. It must be used with DIZET and BHANG. Double-precision variables have been used 
throughout the program in order to obtain maximum accuracy, which is especially important 
for resonance physics. In all, the package {zpl Tj^^, DIZET and BHANG) contains about 11500 

rp 

lines of FORTRAN code. A block diagram of zpl Tj^^^ is shown in fig. |T^. 

The following routines are normally called in the initialization phase of programs using the 

rp 

zpl Tj^^ package. Normally they are called in the order listed below. 



6.1 Subroutine ZUINIT 

Subroutine ZUINIT is used to initialize variables to their default values. This routine must be 

rp 

called before any other zpl Tj^ji routine. 
CALL ZUINIT 



6.2 Subroutine ZUFLAG 

rp 

Subroutine ZUFLAG is used to modify the default values of flags which control various zpl Tj^j^ 
options. 

CALL ZUFLAG (CHFLAG,IVALUE) 
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ZFTEST 



ZUINIT 
1 



ZUFLAG 



DIZET 



I 



ZUWEAK 



ZUINFO 



ZUCUTS 



Standard Model 
interfaces 



ZUTHSM 
ZUTPSM 



Model-independent interfaces 


Effective couplings 


Partial widths 


S-matrix 


ZUXSA 
ZUXSA2 
ZUTAU 


ZUXSEC 


ZUSMAT 



INDF= 11 



BHANG 



ZCUT 



BHAINI 



EWINIT 



BHACOU 



ICUT= 



EWCOUP 



ICUT= 



BHACUT 



-1,0,1 



SCUT 



Figure 14: The structure of zfI T^h- ICUT=-1 gives observables without any cuts, ICUT=0,1 
with cuts. 
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Input Arguments: 

rji 

CHFLAG is the character identifier of a zfI Tj^ji flag (see table P). 
IVALUE is the value of the flag. 



CHFLAG 


IVALUE 


CHFLAG 


IVALUE 


CHFLAG 


IVALUE 


AFBC 


1 


1 


ALPH 








ALST 


1 





AMT4 


3 


3 


BORN 








BOXD 





1 


CONV 





1 


FINR 


1 


1 


F0T2 


1 


2 


GAMS 


1 


1 


INCL 


1 


0;1 


INTF 


1 


1 


PART 








POWR 


1 


1 


PRNT 







QCDC 


1 


2 


QCD3 


1 


1 


VPOL 


3 


3 


WEAK 


1 


1 















Table 9: Flag settings for zpl Tj^^; the values shown are: in the first column the default, recom- 
mended settings optimized for LEP I physics and in the second the 'best' settings, recommended 
for use in a broader energy region. 



Possible combinations of CHFLAG and IVALUE are listed below: 

CHFLAG= ' AFBC ' Controls the calculation of the forward backward asymmetry for interface 
ZUTHSM. 

IVALUE=0 Asymmetry calculation is inhibited (can speed up the program if asymmetries 
are not desired). 

IVALUE=1 (default) Asymmetry calculation is done. 



CHFLAG='ALPH' Controls the calculation of a,(g2, A,^) for i?QCD ( ^251 |3]7|)-(|3]8|) and for 
the 0{aas) corrections (section p.4|) ; see also flag ALST. 

IVALUE=0 (default) ALPHA4 is used to calculate a^, where is calculated according 
to table user input is deflned by ALST. 

IVAHJE=1 ALPHAS relies on 5 quark flavors and is used to calculate a^. 

CHFLAG= ' ALST ' Determines how the strong coupling constant runs as a function of s in 
i?QCD, ( |3.25| ). Form factor corrections 0{aas) are calculated corresponding to ALST=0. 
ALFAS is input by the user in calls to ZUWEAK, ZUTHSM and ZUTPSM. 



IVALUE=0 Calculation of =ALPHAn(g^, ky^) corresponding to (|3.81| ) with = M'z for 



-i\L.rn.i\i.i.\Li uuiitjspuiiumg lu ^|o.oi| ; wiLii "'^^ 

the Z width and with = s for the cross sections. Where n is determined according 

= 4, A^ 

ALFAS. 



to the flag ALPH. For n = 4, Aj^ is calculated from the input A^= 



=ALFAS, while 



for n 



,(5) 



IVALUE=1 (default) In -Rqcd, the strong coupling constant runs as follows: Ois{q^) = 
a,(M|) [ ALPHAn(g2, A^^) / ALPHAn (M|, A^) ], with a,(M|^ -atu-ac w — aH 
is flxed at 185 MeV forn = 4 and at 122 MeV for n = 5. 



ALFAS. Here, A^ 
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CHFLAG= ' AMT4 ' Controls calculation of leading 0{a'^mf) terms for Ar and the weak form 
factors as discussed in section ^]4| 



IVALUE=0 No resummation. Ar, p, k as introduced in sections |3n|j P-2| , |3.3| are calculated 
to order 0{a) , with possible inclusion of the 0{aas) corrections depending on 
another flag; 

IVALUE=1 Leading 0{a^mf) corrections are included in Ap and 6p, ( ^.71[ ), ( ^.7^^! ), while 
the 0{aas) terms and Xq, X are neglected there. The latter remain in the remain- 
der part of the form factors; 

IVALUE=2 Common resummation of leading 0{a'^mf) and 0{aas) terms; only the Xq, X 
are neglected; 

IVALUE=3 (default) Resummation as described in [^2], |T8|, |20[; Ap and 5p as defined in 

CHFLAG= ' BORN ' Controls calculation of QED and Born observables. 

IVALUE=0 (default) QED convoluted observables. 
IVALUE=1 Non-convoluted 'effective' Born observables. 



CHFLAG= ' BOXD ' Determines if the ZZ and WW box contributions (see fig. [T^) are calculated. 

IVALUE=0 (default) No box contributions are calculated. 
IVALUE=1 ZZ and WW box diagrams are calculated. 

CHFLAG= ' CONV ' Controls the energy scale of running a. 

IVALUE=0 (default) For WEAK=1, a is calculated at the energy scale s and for WEAK=0 it 
is not running at all. 

IVAHJE=1 a is calculated at the energy scale s' and convoluted. 
CHFLAG='FINR' Controls the calculation of final-state radiation. 
IVALUE=0 Final-state radiation is included as in ( |1 . 1 1|) . 

IVALUE=1 (default) Include complete treatment of final-state radiation with common 
soft-photon exponentiation as in ( p..lO|) , ( p.l|) . 

CHFLAG= ' F0T2 ' Controls second-order leading log and next-to-leading log QED corrections. 

IVALUE=0 Second-order QED corrections are not included. 

IVAHJE=1 (default) Second-order QED corrections are included as described in |26| ; 
constant terms O(tt^) omitted. 



IVALUE=2 Second-order QED corrections are included as described in |2T|, |26[| . 
CHFLAG= ' GAMS ' Controls the s dependence of Qz, the Z- width function, introduced in 
IVALUE=0 Forces Qz to be constant. Propagator definition (|3.9| ) is used. 



To reproduce results presented in |31 , one has to choose AMT4=-1 together with VPDL=2. This combination 
of flags gives a reasonable result only at LEP I energies. 



42 



IVALUE=1 (default) Allows Qz to vary as a function of s as in 

CHFLAG=' INCL' Influences the treatment of final-state bremsstrahlung exclusively for quarks 
and hadrons. 

IVALUE=0 Same as FINR=1. 

IVALUE=1 (default) For quarks and hadrons, final-state bremsstrahlung is treated as 
with FINR=INTF=0. 

CHFLAG= ' INTF ' Determines if the 0{a) initial-final state QED interference terms are calcu- 
lated. These terms are very small near the Z peak; however, they can become significant 
if severe kinematic cuts are applied. 

IVALUE=0 The interference term is ignored. 
IVALUE=1 (default) The interference term is included. 

CHFLAG= ' PART ' Controls the calculation of various parts of Bhabha scattering. 

IVALUE=0 (default) Calculation of full Bhabha cross section and asymmetry. 
IVALUE=1 Only s channel. 
IVALUE=2 Only t channel. 
IVALUE=3 Only s-t interference. 

CHFLAG='PRNT' Controls ZUWEAK printing. 

IVAHJE=0 (default) Printing by subroutine ZUWEAK is suppressed. 
IVALUE=1 Each call to ZUWEAK produces some output. 

CHFLAG= ' POWR ' Controls inclusion of final-state masses in kinematical factors. 

IVALUE=0 Lepton and light-quark masses are set to zero in ( p.l7|) , (|3.22| ), (|3.55|) , ( p.52|) , 
(U); rric = 1.5, m,, = 4.5 GeV. 

IVALUE=1 (default) Lepton and light-quark masses as taken in the calculation of vacuum 
polarization. In combination with VP0L=1: rriu = .062, = .083, = .215 GeV, 
and with VP0L=2,3: m„ = .04145, = .04146, = .15 GeV. 

CHFLAG= ' QCDC ' Controls how 0{aas) corrections related to the t-quark mass are treated 
within weak form factors. The leading 0{aasmf) term of this QCD correction is expli- 
citly given in ( p.73| ). 

IVALUE=0 0{aas) corrections to weak form factors are not calculated. This setting must 



be used for numerical comparisons with the tables shown in [l3T| , pAj ^0 



IVALUE=1 (default) They are included and determined with a fast approximate calcula- 
tion as described in [^. It should be noted that this approximation is only valid 
for s < m\. 

IVALUE=2 Same as 1 except that exact calculations of the Feynman diagrams are per- 
formed as in No restriction on s. 



An alternative calculation |pa agrees numerically to 12 digits. 
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CHFLAG='QCD3' Controls the inclusion of the in (^7^), (^^ ). 

IVALUE=0 (default) This term is not included. 
IVALUE=1 The calculation is made to O(a^). 

CHFLAG= ' WEAK ' Determines if weak loop calculations are to be performed. 

IVALUE=0 No weak loop corrections to the cross sections are calculated and weak pa- 
rameters are forced to their Born values, i.e. Pe/ = ^e./,e/ = 1- 

IVALUE=1 (default) Weak loop corrections to the cross sections are calculated. 

CHFLAG= ' VPOL ' Controls, which parametrization of the hadronic vacuum polarization contri- 
bution cthad to the photon propagator ( |3.33|) is used. Three different parametrizations 
are available. 



IVAHJE=1 Selects a parametrization taken from [p9| , |61[] . 

IVALUE=2 Quarks are treated like leptons and their effective masses are as in the second 
set quoted in the description of the flag POWR. This choice was used to obtain the 
results presented in [00. 



IVALUE=3 (default) Selects a parametrization that uses the hadronic vacuum polariza- 
tion calculations described in [ITOi. 



6.3 Subroutine ZUWEAK 

Subroutine ZUWEAK is used to perform the weak sector calculations. These are done internally 



with DIZET |10|. The routine calculates a number of important electroweak parameters ( i.e. 
sin^^ty, the partial Z widths, fermionic vacuum polarization. Fa, and weak form factors for 
the cross section), which are stored in common blocks for later use (see appendix A). If any 

zpl Tj^ji flags are to be modified this must be done before calling ZUWEAK. 



CALL ZUWEAK (ZMASS , TMASS , HMASS , ALFAS) 



Input Arguments: 

ZMASS is the Z mass (Mz) in GeV. 

TMASS is the top quark mass (nit) in GeV, [10-400]. 

HMASS is the Higgs mass (M^) in GeV, [10-1000]. 

ALFAS is the value of the strong coupling constant (a^) at = M| (see also flag ALST). 

A tremendous saving in computing time can be realized by performing weak sector calcu- 
lations only once during initialization of the zpl Tj^j^ package. This is possible because weak 
parameters are nearly independent of s near the Z peak, e.g. ~ lns/M|. 

^"In addition, the value of the AMT4 flag was set to —1. 
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6.4 Subroutine ZUCUTS 



Subroutine ZUCUTS is used to define kinematic and geometric cuts for eacli fermion cliannel. 
In terms of the internal structure of zfI Tj^^, this routine is used to select the appropriate 
QED calculational chain. 

CALL ZUCUTS ( INDF , ICUT , ACOL , EMIN , S_PR , ANGO , ANGl ) 

Input Arguments: 

INDF is the fermion index (see table p!0|). 





Final- 


INDF 


state 




fermions 





vv 


1 


e+e^ 


2 




3 




4 


UU 


5 


dd 


6 


cc 


7 


ss 


8 


tt 


9 


bb 


10 


hadrons 


11 


Bhabha 



Table 10: Indices used by zfI Tj^ji interface routines to select the final-state fermion pair. Note 
that INDF=1 returns only s-channel observables, INDF=8 always returns zero, and INDF=10 
indicates a sum over all open quark channels. 

ICUT controls the kinds of cuts {chain) to be used. 

=-1 : (default) no cuts at all are to be used (fastest). 

= : allows for a cut on the acoUinearity of the // pair and the minimum energy of both 
fermion and antifermion. 

= 1 : allows for a cut on the minimum invariant mass of the // pair. 

ACOL is the maximum acoUinearity angle (^™'^^) of the // pair in degrees (ICUT = 0). 

EMIN is the minimum energy (i?™™) of the fermion and antifermion in GeV (ICUT = 0). 

S_PR is minimum allowed invariant // mass {s') in GeV (ICUT = 1). This is related to the 
maximum photon energy by (|2.5|) , (|2.6| ). 

ANGO (default = 0°) is the minimum polar angle in degrees of the final-state antifermion. 

ANGl (default = 180°) is the maximum polar angle (-(9) in degrees of the final-state antifermion. 
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6.5 Subroutine ZUINFO 

Subroutine ZUINFO prints the values of zpl flags and cuts. 
CALL ZUINFO (MODE) 

Input Argument: 

MODE controls the printing of zpl Tj^^ flag and cut values. 

=0 : Prints all flag values. 
=1 : Prints all cut values. 

7 Interface Routines of zpl 

Each calculational branch of zpl Tj^^ has corresponding interfaces. These interfaces will be 
described below. For the Standard Model branch the cross section and asymmetry interface is 
subroutine ZUTHSM, while for the tau polarization it is subroutine ZUTPSM. Subroutines ZUXSA, 
ZUXSA2 and ZUTAU are interfaces for the effective coupling's branch. The interfaces for the 
partial widths and S-matrix branches are ZUXSEC and ZUSMAT, respectively. 

Note that subroutine ZUWEAK must be called prior to any of the interfaces to be described 
below. As a consequence, flags used in this subroutine can influence the calculation of cross 
sections and asymmetries in the interfaces described now. 

7.1 Subroutine ZUTHSM 

Subroutine ZUTHSM is used to calculate Standard Model cross sections and forward-back- 
ward asymmetries as described in section ^. 

CALL ZUTHSM ( INDF , SQRS , ZMASS , TMASS , HMASS , ALFAS , XS* , AFB*) 
Input Arguments: 

INDF is the fermion index (see table p^). 

SQRS is the centre-of-mass energy (y^) in GeV. 

ZMASS is the Z mass {Mz) in GeV. 

TMASS is the top quark mass {rrit) in GeV, [10-400]. 

HMASS is the Higgs mass (M^) in GeV, [10-1000]. 

ALFAS is the value of the strong coupling constant {as) at = M| (see also flag ALST). 

Output Arguments:^ 

XS is the total cross section (ctt) in nb. 

AFB is the forward-backward asymmetry (Afb)- 

^^An asterisk (*) following an argument in a calling sequence is used to denote an output argument. 
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7.2 Subroutine ZUTPSM 

Subroutine ZUTHSM is used to calculate the Standard Model tau polarization and tau polar- 
ization asymmetry as described in section ^. 



CALL ZUTPSM (SQRS , ZMASS , TMASS , HMASS , ALFAS , TAUPOL* , TAUAFB*) 



Input Arguments: 

SQRS is the centre-of-mass energy (y/s) in GeV. 
ZMASS is the Z mass (M^) in GeV. 
TMASS is the top quark mass (mt) in GeV, [40-300]. 
HMASS is the Higgs mass (Mh) in GeV, [10-1000]. 

ALFAS is the value of the strong coupling constant (a^) at = M| (see also flag ALST). 
Output Arguments: 

TAUPOL is the tau polarization (Apoi) of ( |3.67| ). 

TAUAFB is the tau polarization forward-backward asymmetry (A^j^) as defined in ( |3.68 ). 



7.3 Subroutine ZUXSA 



Subroutine ZUXSA is used to calculate the cross section and asymmetry described in section 
as a function of y/s, Mz, Tz, and the weak couplings ([4.8|) , (|4.9|) , (|4.10| ). 



CALL ZUXSA (INDF , SQRS , ZMASS , GAMZ , MODE , GVE , XE , GVF , XF , XS* , AFB*) 



Input Arguments: 

INDF is the fermion index (see table 0), (1:9,11). 

SQRS is the centre-of-mass energy (i/i) in GeV. 

ZMASS is the Z mass (Mz) in GeV. 

GAMZ is the total Z width (Tz) in GeV. 

MODE determines which weak couplings are used: 

=0 : XE (XF) is the effective axial- vector coupling {cja) for electrons (final-state fermions). 

=1 : XE (XF) is the effective weak neutral-current amplitude normalization (p) for elec- 
trons (final-state fermions). 

GVE is the effective vector coupling for electrons (^^). 
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XE is the effective axial- vector coupling (gl) or weak neutral-current amplitude normalization 
(pe) for electrons (see MODE). 

GVF is the effective vector coupling for the final-state fermions {§:[)■ 

XF is the effective axial- vector coupling or weak neutral-current amplitude normalization 
(p/) for the final-state fermions (see MODE). 

Output Arguments: 

XS is the cross section (ctt) in nb. 

AFB is the forward-backward asymmetry {Aps)- 

7.4 Subroutine ZUXSA2 

Subroutine ZUXSA2 is used to calculate the lepton cross section and asymmetry as a function 
of ^/s, Mz, T Zi and the weak couplings assuming lepton universality. This routine is similar 
to ZUXSA except that the couplings are squared. 

CALL ZUXSA2 ( INDF , SQRS , ZMASS , GAMZ , MODE , GV2 , X2 , XS* , AFB* ) 
Input Arguments: 

INDF is the fermion index (see table 0) (1-3,11). 

SQRS is the centre-of-mass energy (i/i) in GeV. 

ZMASS is the Z mass {Mz) in GeV. 

GAMZ is the total Z width {T z) in GeV. 

MODE determines which weak couplings are used: 

=0 : X2 is the square of the effective axial- vector coupling (^^) for leptons. 

=1 : X2 is the square of the effective neutral-current amplitude normalization (p/) for 
leptons. 

GV2 is the square of the effective vector coupling (^^) for leptons. 

X2 is the square of the effective axial-vector coupling (^^) or neutral-current amplitude nor- 
malization (p;) for leptons (see MODE). 

Output Arguments: 

XS is the cross section (ctt) in nb. 

AFB is the forward-backward asymmetry {Afb)- 
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7.5 Subroutine ZUTAU 

Subroutine ZUTAU is used to calculate the r"*" polarization as a function of ^/s, Mz, Tz, and 
the weak couplings (see discussion in section 



CALL ZUTAU (SQRS , ZMASS , GAMZ , MODE , GVE , XE , GVF , XF , TAUPOL* , TAUAFB* ) 



Input Arguments: 

SQRS is the centre-of-mass energy (i/i) in GeV. 

ZMASS is the Z mass (Mz) in GeV. 

GAMZ is the total Z width (Tz) in GeV. 

MODE determines which weak couplings are used: 

=0 : XE (XF) is the effective axial- vector coupling (cja) for electrons (final-state fermions). 

=1 : XE (XF) is the effective weak neutral-current amplitude normalization (p) for elec- 
trons (final-state fermions). 

GVE is the effective vector coupling for electrons (^^). 

XE is the effective axial- vector coupling {g^) or weak neutral-current amplitude normalization 
(pe) for electrons (see MODE). 

GVF is the effective vector coupling for the final-state fermions {§{, )■ 

XF is the effective axial- vector coupling {g[) or weak neutral-current amplitude normalization 
[pf) for the final-state fermions (see MODE). 

Output Arguments: 

TAUPOL is the tau polarization (A^) defined in (|3.67|) . 

TAUAFB is the forward-backward asymmetry for polarized tau's (A^j^) as defined in (|3.68|) . 



7.6 Subroutine ZUXSEC 

Subroutine ZUXSEC is used to calculate the cross section as a function of i/i, Mz, Tz, Fg and 
F/ as was described in section | 



CALL ZUXSEC ( INDF , SQRS , ZMASS , GAMZO , GAMEE , GAMFF , XS*) 



Input Arguments: 

INDF is the fermion index (see table pil|). 

SQRS is the centre-of-mass energy (y^) in GeV. 

ZMASS is the Z mass (Mz) in GeV. 
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GAMZO is the total Z width (F^) in GeV. 
GAMEE is the partial Z decay width (Fe) in GeV. 

GAMFF is the partial Z decay width (F/) in GeV; if INDF=10, GAMFF=r/,. 

Output Argument: 

XS is the cross section (crT)in nb. 

7.7 Subroutine ZUSMAT 

Subroutine ZUSMAT is used to calculate the cross section from the S-matrix approach (see 
section |0) . 



CALL ZUSMAT ( INDF , SQRS , ZMASS , GAMZ , RR , RI , RO , Rl , R2 , RG , XS * ) 



Input Arguments: 

INDF is the fermion index (see table 0), [2,10]. 
SQRS is the centre-of-mass energy (y^) in GeV. 
ZMASS is the Z mass [Mz) in GeV. 
GAMZ is the total Z width (F^) in GeV. 

RR-R2 six parameters in S-matrix approach, (RR, RI, RG, RO, . . .) = {R, I ^r^^r^, . . .), intro- 
duced in i W^ - 



Output Argument: 

XS is the cross section (ctt) in nb. 

Note that the default Z mass and width definitions correspond to (|3^) and thus differ 
from those of the other zpl Tj^^ interfaces. 

8 Comparisons 

In this section, we compare the predictions of the Standard Model branch of zf"^^ with 
other programs. For this comparison we use the following parameter values, unless explicitly 
stated otherwise: Mz = 91.18, rut = 150, Mh = 100 GeV, and as = 0.12. The section is 
broken up into two parts: 

1 . A comparison of the weak mixing angle (with its various definitions) as well as the partial 
and total Z widths. 

2. A comparison of total cross sections and forward-backward asymmetries. 
In the past many comparisons of this sort have been made. In particular an earlier version of 



zpl TjTi^, the ZBIZON [|7T| package, was used in the 1989 Workshop on Z Physics at LEP 1 |^ for 



numerous comparisons |T^, 0, ^ . In addition to these comparisons, others did not explicitly 
include the ZBIZON code |]72|. At that time, the predictions from all of these programs agreed 



to within 0.5%. Since then, new codes have been developed and the quality of several of the 
existing programs has been improved; among these is the zpl Tj^^ package. 
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8.1 Weak Mixing Angles and Partial Z Widths 



Throughout the Standard Model branch of zf^ T^^^ we use the on-shell definition ( |3.15D of the 



weak mixing angle sin Qt^. For the sake of this comparison, we take into account two additional 



definitions of the weak mixing angle: the 'effective' weak mixing angle, s^^, introduced in 



( p.20| ), and the weak mixing angle of the MS renormahzation scheme given below: 



sin^ Qy/ 



2,eff 



Mr 



sin^ 9w, 



(8.1) 



(8.2) 



sin^^^^s 



1 + . o ^ op 



sin^ e 



w 



sin^ e 



w, 



(8.3) 



where 6p has been introduced in ( p. 73 ). In table 11, we compare predictions from zpI^Tj^^ 
with those obtained by W. Hollik ||T^ and G. Degrassi, S. Fanchiotti, A. Sirlin |T8 |. 



The agreement of the different calculations for these three cases is impressive. From the 
table it is apparent that the mixing angles Sy^^ and sin^ 9^ depend to a lesser extent on the 
unknown top and Higgs masses than does sin^^vi^- For a detailed discussion of the different 
approaches see for instance |]73|, ^ ^ . 

When using data to determine an effective weak mixing angle, one must be careful, since 
measurements of mixing angles from different observables may yield results that cannot be 
directly compared. This delicate point was addressed in sections ^]2| and ^TT] [see also ( |3.43| )- 



( p.45|) ] and has been discussed in detail in ||4^. It has been demonstrated that a proper 
formulation of the hard-scattering subprocess and a correct unfolding of the leptonic forward- 
backward, b-quark forward-backward, and tau polarization asymmetries lead to results which 



are very close to each other and to the value of expected from Fg, as defined in ( p. 20 ). 



In table |T^, we compare partial, hadronic and total Z widths with numbers of other authors: 
as W. Hollik |T^ and G. Degrassi, A. Sirlin Shown in the second and third lines are the 



digits which differ from zpl Tj^^ ^ as one can see the deviations are very small. 
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rrit 


Mh 


sin^ Ov/ 


2, off 




100 


100 


0.23056 


0.23362 


0.23351 






44 


52 


7 






62 


66 


42 




500 


0.23266 


0.23447 


0.23438 






53 


37 


43 






73 


51 


29 




1000 


0.23371 


0.23485 


0.23477 






60 


77 


83 






81 


91 





150 


100 


0.22483 


0.23217 


0.23213 






74 


07 


25 






87 


23 







500 


0.22690 


0.23300 


0.23299 






81 


291 


309 






5 


6 


5 




1000 


0.22794 

88 


0.23337 



0.23337 
49 






802 


45 


5 


200 


100 


0.21782 


0.23025 


0.23024 






78 


17 


41 






7 


35 


8 




500 


0.21985 


0.23106 


0.23108 






3 


099 


24 






91 


16 


10 




1000 


0.22088 


0.23142 


0.23144 






7 


37 


63 






96 


54 


9 


250 


100 


0.20919 


0.22786 


0.22785 






21 


1 


808 






25 


800 


96 




500 


0.21118 


0.22865 


0.22866 






22 





89 






26 


79 


76 




1000 


0.21217 


0.22899 


0.22901 






25 


7 


26 






28 


916 


13 



Table 11: Comparison ofsn^Ow, Sy^^ and sin^ 6'Jy^(M|) from zpJ^Tj^^ (first line), Hollik / (I^/ 
(second line), and Degrassi, Fanchiotti, Sirlin J7^ (third line); with flags AMT4=3, QCDC=0, 
QCD3=0, Mz=91.170 GeV, rut and Mh m GeV. 
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nit 


Mh 


















Tb 


Thad 


Ttot 


100 


100 


166.3 


83.42 

4 


83.42 

4 


83.23 


296.0 


382.2 
1 


295.6 


382.2 


377.5 
6 


1733.5 


2482.4 




500 


166.1 


83.29 
30 
3 


83.29 
30 


83.10 


295.3 
2 
2 


381.4 
2 
3 


294.8 


381.4 


376.7 


1729.5 


2477.5 




1000 


166.0 


83.22 
3 

Z 


83.22 
3 


83.03 


294.8 
7 


380.9 
7 

Q 



294.4 


380.9 


376.2 
3 


1727.3 


2474.8 


150 


100 


166.9 
7.0 


83.81 
3 

8 


83.81 
3 


83.62 


298.0 
7.9 
7.9 


384.4 


297.5 


384.4 


376.5 
6 


1740.8 


2492.9 




500 


166.8 


83.68 
9 


83.68 
9 


83.49 


297.2 


383.6 
5 


296.8 


383.6 


375.7 

8 


1736.9 


2488.0 




1000 


166.6 


83.61 

D 


83.61 


83.42 


296.8 
7 

( 


383.1 



296.4 


383.1 


375.3 


1734.7 


2485.3 


200 


100 


167.9 

8 


84.37 

8 


84.37 

8 


84.18 


300.6 

5 


387.4 


300.2 


387.4 


375.3 


1750.9 


2507.4 




500 


167.7 


84.24 

O 


84.24 


84.04 


299.9 

8 

Q 
O 


386.6 


299.4 


386.6 


374.5 


1747.1 


2502.6 




1000 


167.6 

5 


84.16 

2* 


84.16 


83.97 


299.5 

4 
4 


386.2 
1 
1 


299.0 


386.2 


374.1 


1745.0 


2499.9 


250 


100 


169.0 
1 


85.10 
1 

1* 


85.10 
1 


84.91 


304.0 
3.8 


391.4 
3 


303.6 


391.4 


373.7 

5 


1764.0 


2526.3 




500 


168.9 

8 


84.96 
7 

5.0* 


84.96 
7 


84.77 


303.3 
1 


390.6 

5 


302.9 


390.6 


373.0 
2.8 


1760.2 


2521.5 




1000 


168.8 
7 


84.89 

9* 


84.89 


84.69 


302.9 
7 

8 


390.1 


302.5 


390.1 


372.6 

4 


1758.2 


2518.9 



Table 12: Partial and total widths of the Z boson in MeV from zfI Tj^^ (first line), Hollik ^T^J 
(second line) and Degrassi, Sirlin f^D^ (third line). Shown are only the digits which differ, a 
dash means no entry, an asterisk no digit available. Flags as in table\T^, Mz=91.170 GeV, 
rrit and Mh in GeV. 
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8.2 Cross Sections and Asymmetries 

In this section the cross sections and asymmetries for processes ( |1 . ID and ( |1.2[ ) predicted by 



zfI^^Er compared with those of ZSHAPE 2.0 |74[ and of ALIBABA 2.0 



Earher 





751 




i in 



For (|1 . 1| ) , all packages include weak corrections of at least 0{a). In ZSHAPE, QED con- 
tributions of 0(a;^) to the initial-state are calculated exactly, while in the other programs a 
leading-log approximation is used. In addition, all programs include final-state radiation cor- 
rections to 0{a) and common exponentiation of initial- and final-state soft-photon emission. 
Initial-final interference is contained only in ALIBABA and zfI Tj^ji- The additional t-channel 
terms (including s-t interference) , which are necessary in order to calculate Bhabha scattering 

rri 

( |1.2|) , are available in both ALIBABA and zpl Tj^^. In the latter, this is done via the BHANG 

package, which only contains some of the higher-order t-channel QED corrections that have 

T 

been implemented in ALIBABA. On the other hand, zpl Tj^^ contains higher-order weak and 



QCD corrections as explained in section |3.4| , which are not available in the other two pack- 
ages. As can be seen in more detail from the references, the three codes allow for different 
applications of kinematic cuts due to their different theoretical basis. 

Since the other two programs (ALIBABA and ZSHAPE) perform only Standard Model cal- 
culations, we have restricted these comparisons to the corresponding branch of zpl Tj^^. In 
this context, we have used the 'recommended' zfI^Tjt^ji flags of table ||. In addition, we per- 



form comparisons using flag settings shown in table 13, which have been chosen such that the 



corrections realized in zpl Tj^^ most closely resemble that of the other programs. 



CHFLAG 


IVALUE 


CHFLAG 


IVALUE 


CHFLAG 


IVALUE 


AFBC 


1 


1 


ALPH 


1 


1 


ALST 


1 


1 


AMT4 








BORN 








BOXD 





1 


CONV 


1 


1 


FINR 


1 


1 


F0T2 


2 


2 


GAMS 


1 


1 


INCL 








INTF 





1 


PART 








POWR 








PENT 






QCDC 








QCD3 








VPOL 


2 


2 


WEAK 


1 


1 















Table 13: Flag settings in zpl Tjt^r for comparisons of cross sections and asymmetries; the 
values shown are: first column - best agreement with ZSHAPE, second column - best agreement 
with ALIBABA. 

We performed three series of comparisons: 

I. A comparison of the total zpl Tj^^ and ZSHAPE cross sections with an s' cut and no 
angular acceptance cut. 

II. Comparisons of zpl Tj^^ and ALIBABA muon pair production cross section and forward- 
backward asymmetry with i?™™, and angular acceptance cuts. 



HI. As above, except that here the comparison is done for Bhabha scattering. 
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For case I, we have compared quark and muon cross sections over a large energy range, 
[10-100] GeV. Due to ZSHAPE limitations, only cross-section comparisons up to LEP I energies 



can be performed. In fig. |T^, we show the ratio of hadron and muon cross sections as a 



function of the centre-of-mass energy for different values of s'. As can be seen in the figure, 
the agreement at LEP I energies is excellent even though some higher-order weak corrections 
are not realized in ZSHAPEQ. As the energy decreases from LEP I, deviations begin to appear, 
which reach 2% in magnitude for the hadronic cross section. For muons, the deviation goes 
to 1% when s' is at the kinematic limit and 0.5% otherwise, approaching in the latter case 
nearly exact agreement at small energies where pure QED dominates. 

A considerable improvement in the agreement of these two programs can be realized 
through a judicious choice of zpl Tj^^ flags (see table |1^). If various enhancements to zpl Tjr^j^, 
which have been realized since the 1989 workshop, are inhibited, then a dramatic improvement 
in the agreement of the two programs is observed. These enhancements are mainly concerned 
with QCD corrections, higher order weak corrections and the handling of light-quark thresh- 
olds. As can be seen from fig. [16], the disagreement shrinks to 0.1% for both muon and hadron 
production cross sections and for different cuts. 

In case II, we compare the predictions of zpl Tj^^ with ALIBABA for the muon produc- 
tion cross section and forward-backward asymmetry with cuts on minimum fermion energy 
[E™™) and acoUinearity {S,^^^). We restrict this comparison to muons since ALIBABA has no 
hadron option. In figs. |1^ and |T8|, we contrast the predictions for these programs using the 
'recomended' zpi^Tj^p^^ flags and another set of flags (table |TB[) chosen to minimize the differ- 
ences in the calculations performed by these programs. The differences in the predictions, for 
the value of mt chosen, is minor in both cases. For the 'recommended' flags, over the large 
energy interval covered in fig. |1^, the cross sections agree to within 0.7%; at LEP I energies the 
agreement is good to within 0.2%. For the asymmetry, the difference of the two predictions is 
smaller than 0.2% over the full energy range and well within 0.1% at LEP I energies. 

For case III, we compare in fig. |l^ the cross-section ratio and forward-backward asymmetry 
for Bhabha scattering for zpl Tjr^j^ (via BHANG) and ALIBABA with the same cuts as described 
above for case II. Since BHANG contains several approximations adapted to applications at 
LEP I, we restrict the energy range of the comparison correspondingly. As may be seen from 
the figure, the programs agree to within 1.5% for the cross-section ratio and within 1% for the 
asymmetry difference. 



rp 

corrections connected with the t quark. The agreement of the program with zpl Tj^j^^ with the corresponding 



Very recently, a new version of ZSHAPE was developed, which now also contains higher-order electroweak 
ections connected with the t quark. The ag 
flag settings, is not worse than shown here 
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Figure 15: Ratio of cross section predictions from zfI T^ji and ZSHAPE for muon and hadron 
production, as a function of the centre- of -mass energy, for three different values of A — 1 — 
s^in/s; both programs with their 'recommended' choice of flags. 
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Figure 16: Ratio of cross-section predictions from zpl Tj^^^ and ZSHAPE for muon and quark 
production, as a function of the centre- of -mass energy, as in fig. [I^, hut here flags are chosen 
such that the theoretical assumptions of both programs are as similar as possible. 
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Figure 17: Ratio of cross sections, , and difference of forward-backward asymmetries, Ap^, 
as predicted by zpl Tj^^ anc? ALIBABA, as a function of the centre- of -mass energy. An acceptance 
cut of 44° < < ISff' , an acollinearity cut of ^ < 25° and a muon energy cut of E 

5 GeV have been employed. The comparison is made for both the 'recommended' zfI Tj^j^ f^cLQ 
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Figure 18: Same as in fig. [7^, for LEP I energies explicitly. 
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8.3 Conclusions 

There exists a wealth of programs, of the semi-analytical and Monte Carlo variety, which can 
make predictions for the fermion pair production process in e~^e~ collisions. The programs 
are varied in both their theoretical accuracy ( i.e. the order to which the calculations are 
performed) and in the cuts that one may apply. zfI Tj^^ is a semi-analytic program with 

large inherent flexibility in both these respects. With a judicious selection of zpl Tj^^ flags, 
agreement with other programs to the level of 0.5% and better have been reached with the 
exception of Bhabha scattering where the agreement is slightly worse. 
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Note added in proof: 

Recently, a new package has been published, which allows to calculate the one-loop electroweak 
radiative corrections to two- fermion production near the Z resonance |8l|. The test run output for 
the unpolarized muon production cross section shown in Sample 2 (p. 62 of ||8l|) reproduces the 
effective Born approximation of zpI^T^^ within 0.01% at the resonance, and within 0.02% in the full 
energy range covered by Sample 2. The corresponding b-quark production cross section shown in 

rp 

Sample 3 does not agree with the zpl Tj^ji results. 
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A Common Blocks 

A.l zpl T^ji Common Blocks 

zpl Tj^^ common blocks of potential interest to the user are documented here. 

COMMON /ZUPARS/QDF,QCDC0R,QCDC0B,ALPHST,SIN2TW,S2TEFF(0:11) , 
& WIDTHS (0:11) 

The common block /ZUPARS/ contains some zpI^Tj^ff parameters: 
QDF is the final-state radiation factor |^ introduced in (|1 . 1 1| ) and ( p.23| ). 
QCDCOR is a QCD correction for all final quark states except hb defined in (|3.79| ). 
QCDCOB is a QCD correction for hh final states defined in ( |3.8CI|) . 
ALPHST is as{Mz) and is calculated as defined by fiag ALST. 
SIN2TW is sin^6'vF as in (gTT|). 

S2TEFF are the values of for each fermion channel (see (|3.21|) and table |iy). Note that 
S2TEFF(10:11) are not defined. 



WIDTHS are the partial decay widths ( |3.17| ) of the Z for fermion channels defined in table |T0 
(WIDTHS (11) is the total Z width). 

C0MM0N/EWF0RM/XALLCH(5,4) ,XFOTF 
COMPLEX* 16 XALLCH,XFOTF 

Electroweak form factors are stored in the common block /EWFORM/. 

XALLCH(I,J) contains the form factors p, Kg, and Kg/ (J = 1-4) for neutrinos, leptons, 
u and c quarks, d and s quarks, and b quarks (I = 1-5). These have been introduced in 

(EH)- 

XFOTF is 1 + Aa(s) as used in (|03|) . 

COMMON/ZFCHMS/ALLCH (0 : 1 1 ) , ALLMS (0:11) 

The common block /ZFCHMS/ contains the charges and masses of the fermions (see 
table 0). 

ALLCH the fermion charges. 
ALLMS the fermion masses. 

Note that ALLCH (10) and ALLMS (10) are undefined. 

We would also like to mention that the variables FAA, FZA, FZZ which are introduced 
as DATA in subroutine EWCOUP allow us to switch on/off the 77, 7^, ZZ parts of the cross 
sections, respectively. 
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A. 2 DIZET Common Blocks 

Two DIZET internal common blocks of potential interest to the user are documented here. 
COMMON/CDZRKZ/ARROFZ (0 : 10) , ARKAFZ (0 : 10) , ARVEFZ (0:10), ARSEFZ (0 : 10) 



Weak form factors, and kJ, and vector couplings, for partial Z widths ( p.l7|) as 
calculated in subroutine ZWRATE. The indices correspond to those of table llOl 



ARROFZ pj, introduced in (|3.17|) 
ARKAFZ introduced in (|3T^ ) 



ARVEFZ as defined in ( pTTOl) . 



ARSEFZ Effective weak mixing angles SyJ as in ( |3.21|) . 

Note that the 10th element of these arrays is undefined. 

COMMON/CDZXKF/XROKF 

This variable is the ratio of two different definitions of the weak mixing angle as defined 
in (pD: 

3?e XROKF = 

sin Ow 



B Subroutine ZFTEST 

The zpi^Tj^^ distribution package includes subroutine ZFTEST which serves essentially three 
purposes: 

1. It is an example of how to use zpl Tj^^. 

2. It is an internal consistency check of the different zpl Tj^^ branches. 

3. It allows one to check that zpl Tj^n has been properly installed on the machine. 

The routine creates a table of cross sections and asymmetries as a function of ^/s near the Z 
peak. 

To run ZFTEST the user needs to create the following main program: 

PROGRAM ZFMAIN 
CALL ZFTEST 
END 

After compiling and linking it to ZFITR4_5, DIZET and BHANG the results presented in 
appendix p.2| should be obtained. The corresponding Fortran files may be found at ZFIT- 
TEROCERNVM. 
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B.l Subroutine ZFTEST 



SUBROUTINE ZFTEST 

* ========== ====== 

* 

* SUBR. ZFTEST 
* 

* Example program to demonstrate the use of the ZFITTER package. 
* 

* 

IMPLICIT REAL*8(A-H,0-Z) 

COMPLEX* 16 XVPOL 

DIMENSION XS(0:11,5) ,AFB(0:11,4) ,TAUP0L(2) ,TAUAFB(2) 

* 

* constants 
* 

PARAMETER (GMU=1 . 166388D-5 , ALFAI=137 . 0359895D0 , ALFA=1 . DO/ALFAI , 
+ C0NS=1.D0) 
PARAMETER (ZMASS=91 . 175D0 , TMASS=140 . DO , HMASS=300 . DO) 
PARAMETER (AME=0 . 51 lD-3 , AMU=0 . 106D0 , AMT=1 . 784D0 , ALFAS= . 120D0) 
PARAMETER (RSMN=87 . DO , DRS=1 . DO , NRS=9) 
PARAMETER(ANG0=35D0,ANG1=145D0) 

PARAMETER(QE=-1.D0,AE=-.5D0,QU= 2.D0/3.D0,AU= .5D0, 
+ QD=-1.D0/3.D0,AD=-.5D0) 

* 

* ZFITTER common blocks 
* 

COMMON /ZUPARS/QDF,QCDC0R,QCDC0B,ALPHST,SIN2TW,S2TEFF(0:11) , 

& WIDTHS(0:11) 

COMMON /CDZRKZ/ARROFZ (0:10), ARKAFZ (0 : 10) , ARVEFZ (0:10), ARSEFZ (0:10) 
COMMON /EWF0RM/XALLCH(5,4) ,XFOTF 
COMPLEX* 16 XALLCH.XFOTF 

* 

* 

* 

* initialize 
* 

CALL ZUINIT 

* 

* set ZFITTER flags and print flag values 
* 

CALL ZUFLAG ( ' PRNT ' , 1 ) 
CALL ZUINFO(O) 

* 

* do weak sector calculations 
* 

CALL ZUWEAK (ZMASS , TMASS , HMASS , ALFAS) 

* 
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* define cuts for fermion channels and print cut values 
* 

CALL ZUCUTSC 1 ,0, i5.D0,20.D0,0.D0,ANG0,ANGl) 
CALL ZUCUTSC 2 , , 15 . DO , 20 . DO , . DO , ANGO , ANGl) 
CALL ZUCUTSC 3 , , 15 .DO , 20 .DO , .DO , ANGO , ANGl) 
CALL ZUCUTSCll, 0,15. DO, 20.D0,0.D0, ANGO, ANGl) 
CALL ZUINFOCD 

* 

* maie table of cross sections and asymmetries 
* 

PI = DACOSC-l.DO) 
GAMZ = WIDTHS Cll)/1000. 
GAME = WIDTHS C 1)/1000. 
GAE = SQRTCARR0FZCl))/2. 
GVE = ARVEFZC1)*GAE 
DO I = 1,NRS 

RS = RSMN+REALCI-1)*DRS 

* table header 

PRINT SQRTCS) = ' ,REALCRS) 

PRINT *,' < Cross Section 

+ ' < Asymmetry >',' < TauPol >' 

PRINT *,'INDF ZUTHSM ZUXSEC ZUXSA ZUXSA2 ZUSMAT' , 
+ ' ZUTHSM ZUXSA ZUXSA2 ZUTPSM ZUTAU' 

* loop over fermion indices 

DO INDF = 0,11 
S=RS**2 

* standard model interf . ClNTRF=l) 

CALL ZUTHSM C INDF , RS , ZMASS , TMASS , HMASS , ALFAS , 
+ XSCINDF,1) ,AFBCINDF,1)) 

IFCINDF.EQ.3) CALL ZUTPSM CRS, ZMASS, TMASS, HMASS, ALFAS, 
+ TAUPOL CD, TAUAFB CD) 

* cross section interf. ClNTRF=2) 

GAME = WIDTHS CINDF)/1000. 
IFClNDF.EQ.ll) GAME = WIDTHSC D/IOOO. 

CALL ZUXSEC C INDF , RS , ZMASS , GAMZ , GAME , GAME , XS C INDF , 2 ) ) 

* cross section & forward — backward asymmetry interf. ClNTRF=3) 

IFCINDF.NE.O .AND. INDF. NE. 10) THEN 
GAE = SQRTCARR0FZClNDF))/2. 
GVE = ARVEFZCINDF)*GAF 
IFClNDF.EQ.il) THEN 

GAF = SQRTCARR0FZCl))/2. 

GVF = ARVEFZC1)*GAF 
ENDIF 

CALL ZUXSA C INDF , RS , ZMASS , GAMZ , , GVE , GAE , GVF , GAF , 
+ XSCINDF,3) ,AFBCINDF,3)) 

ENDIF 

* tau polarization interf. ClNTRF=3) 

IFCINDF.EQ.3) CALL ZUTAU CRS , ZMASS , GAMZ , , GVE , GAE , GVF , GAF , 

+ TAUPOL C 2 ) , TAUAFB C 2 ) ) 

* cross section & forward — backward asymmetry interf. for gv**2 and 
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* ga**2 (IBRA=4) 

IF((INDF.GE.l .AND. INDF.LE.3) .OR. INDF.EQ.ll) THEN 
GVF2 = GVF**2 
GAF2 = GAF**2 

CALL ZUXSA2 (INDF , RS , ZMASS , GAMZ , , GVF2 , GAF2 , 
+ XS(INDF,4) ,AFB(INDF,4)) 

ENDIF 

* S-matrix interf . (INTRF=5) 

* Parameters are fitted with code FITSMA FORTRAN. 

IF (INDF. EQ. 2 .OR. INDF. EQ. 10) THEN 
IF(INDF.EQ.2) THEN 
AMZS =91.1413200 
GAMZS= 2.48354D0 
RR=0.14159D0 
RI=0.15092D0 
RG=1.13310D0 
ELSE 

AMZS =91.1412600 
GAMZS= 2.4840600 
RR=2 . 9402800 
RI=3. 1512500 
RG=2.81462D0 
ENOIF 

C The parameters AMZS, GAMZS, RR, RI, RG correspond to: 



C AMZS = ZMASS-GAMZ**2/2.00/ZMASS 

C GAMZS = GAMZ -GAMZ**3/2 .00/ZMASS**2 

C VE = -. 500+2. 00*SIN2TW 

C VU = 0.500-4. 00/3. 00*SIN2TW 

C VO = -0.500+2. 00/3. 00*SIN2TW 

C AKAPPA = GMU*AMZSl*AMZSl/(SqRT(2.00)*2.00*PI*ALFA) 

C XVPOL = l.O0/(2.O0-XF0TF) 

C IF(IN0F.LE.3) THEN 

C RZ = CONS*AKAPPA**2*(AE**2+VE**2)**2*(1.00+.7500*ALFA/PI) 

C SZ = C0NS*AKAPPA*VE**2*(l.O0+.75O0*ALFA/PI) 

C RG = C0NS*COABS(XVP0L)**2*(l.O0+.75O0*ALFA/PI) 

C ELSE 

C RZ = C0NS*AKAPPA**2*(AE**2+VE**2)*3.O0* 

C + (2.00*(AU**2+VU**2)+3.00*(AO**2+VO**2)) 

C SZ = C0NS*AKAPPA*3.O0*QE*VE*(2.O0*VU*QU+3.O0*VO*QO) 

C RG = C0NS*COABS (XVPOL) **2*3. 00*11. 00/9. 00 

C ENOIF 

C RR ~ RZ 

C RI = RZ+2.00*SZ*OREAL (XVPOL) 
c 



RI = OdO 
R2 = OdO 
R3 = OdO 

CALL ZUSMAT ( INOF , RS , AMZS , GAMZS , RR , RI , RO , RI , R2 , RG , XS ( INOF , 5 ) ) 
ENOIF 

* results 
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IF(INDF.Eq.O) THEN 

PRINT 9000, INDF, (XS(INDF, J) , J=l ,2) 
ELSEIFdNDF.EQ.l .OR. INDF. EQ. 11) THEN 

PRINT 9010, INDF, (XS (INDF, J) , J=l ,4) ,AFB(INDF,1) , 
+ (AFBCINDF, J) , J=3,4) 

ELSEIF(INDF.EQ.2) THEN 

PRINT 9005, INDF, (XS (INDF, J) , J=l ,5) ,AFB(INDF,1) , 
+ (AFBdNDF, J) , J=3,4) 

ELSEIF(INDF.Eq.3) THEN 

PRINT 9015, INDF, (XS(INDF, J) , J=l ,4) ,AFB(INDF,1) , 
+ (AFBCINDF, J) ,J=3, 4) , (TAUPOL(J) ,J=1,2) 

ELSEIF(INDF.EQ.IO) THEN 

PRINT 9025, INDF, (XS (INDF , J) , J=l ,2) ,XS(INDF,5) 
ELSE 

PRINT 9020, INDF, (XS(INDF, J) , J=l ,3) ,AFB(INDF,1) , 
+ AFBdNDF, 3) 

ENDIF 
ENDDO 
PRINT * 
ENDDO 
RETURN 
9000 FORMAT (IX, 14, 2F8. 4) 
9005 FORMAT (IX, 14, 9F8. 4) 
9010 FORMAT (IX , 14 , 4F8 . 4 , 8X , 3F8 . 4) 
9015 FORMAT ( IX , 14 , 4F8 . 4 , 8X , 5F8 . 4) 
9020 FORMAT ( IX , 14 , 3F8 . 4 , 16X , 2F8 . 4) 
9025 F0RMAT(1X,I4,2F8.4,16X,F8.4) 
* END ZFTEST 

END 
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B.2 ZFTEST Results 



** 


This 


is ZFITTER version 4.5 


** 


** 




92/04/19 


** 




** 


The authors of the ZFITTER package are: 


** 


** 






** 


** 


D . Bardin 


(Dubna) 


** 


** 


M.Bilenky 


(Dubna) 


** 


** 


A . Chizhov 


(Dubna) 


** 


** 


A . Olshevsky 


(Dubna) 


** 


** 


S .Riemann 


(Zeuthen) 


** 


** 


T.Riemann 


(Zeuthen) 


** 


** 


M.Sachwitz 


(Zeuthen) 


** 


** 


A . Sazonov 


(Dubna) 


** 


** 


Yu.Sedykh 


(Dubna) 


** 


** 


I . Sheer 


(UC San Diego) 


** 


** 






** 



** Questions and comments to ZFITTER@CERNVM.CERN.CH ** 
ZUINIT> ZFITTER defaults: 



ZFITTER flag values: 



AFBC: 


1 


ALPH: 





ALST: 


1 


AMT4: 


3 


BORN: 





BOXD: 





CONV: 





FINR: 


1 


F0T2: 


1 


GAMS: 


1 


INCL: 


1 


INTF: 


1 


DUMY: 





PART: 





POWR: 


1 


PRNT: 





qCD3: 


1 


QCDC: 


1 


VPOL: 


3 


WEAK: 


1 



ZFITTER cut values: 



IMDF 


I CUT 


ACOL 




EMIN 




S_PR 


ANGO 


ANGl 





-1 





00 





0000 





0000 





00 


180.00 


1 


-1 





00 





0000 





0000 





00 


180.00 


2 


-1 





00 





0000 





0000 





00 


180.00 


3 


-1 





00 





0000 





0000 





00 


180.00 


4 


-1 





00 





0000 





0000 





00 


180.00 


5 


-1 





00 





0000 





0000 





00 


180.00 


6 


-1 





00 





0000 





0000 





00 


180.00 


7 


-1 





00 





0000 





0000 





00 


180.00 


8 


-1 





00 





0000 





0000 





00 


180.00 


9 


-1 





00 





0000 





0000 





00 


180.00 


10 


-1 





00 





0000 





0000 





00 


180.00 


11 


-1 





00 





0000 





0000 





00 


180.00 



ZFITTER flag values: 



AFBC: 


1 


ALPH: 





ALST: 


1 


AMT4: 


3 


BORN: 





BOXD: 





CONV: 





FINR: 


1 


F0T2: 


1 


GAMS: 


1 


INCL: 


1 


INTF: 


1 


DUMY: 





PART: 





POWR: 


1 


PRNT: 


1 


QCD3: 


1 


QCDC: 


1 


VPOL: 


3 


WEAK: 


1 
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ZMASS = 
HMASS = 
ALPHST = 
QCDCOR = 



91.17500 
300.00000 
0.12000 
1.03954 



TMASS = 
ALFAS = 

SIN2TW = 
QCDCOB = 



140.00000 
0.12000 
0.22817 
1.04020 



CHANNEL 



WIDTH 



nu , nubar 


166 


6 


e+,e- 


83 


6 


mu+ , mu- 


83 


6 


tau+,tau- 


83 


4 


u , ubar 


296 


6 


d , dbar 


382 


9 


c , cbar 


296 


2 


s , sbar 


382 


9 


t ,tbar 








b , bbar 


375 


7 


hadron 


1734 


2 


total 


2484 


7 



ZFITTER cut values: 



DP 


I CUT 


ACOL 




EMIN 




S_PR 


ANGO 


ANGl 





-1 





00 





0000 





0000 





00 


180 


00 


1 





15 


00 


20 


0000 





0000 


35 


00 


145 


00 


2 





15 


00 


20 


0000 





0000 


35 


00 


145 


00 


3 





15 


00 


20 


0000 





0000 


35 


00 


145 


00 


4 


-1 





00 





0000 





0000 





00 


180 


00 


5 


-1 





00 





0000 





0000 





00 


180 


00 


6 


-1 





00 





0000 





0000 





00 


180 


00 


7 


-1 





00 





0000 





0000 





00 


180 


00 


8 


-1 





00 





0000 





0000 





00 


180 


00 


9 


-1 





00 





0000 





0000 





00 


180 


00 


10 


-1 





00 





0000 





0000 





00 


180 


00 


11 





15 


00 


20 


0000 





0000 


35 


00 


145 


00 



72 



SQRT(S) = 87.0000000 





<- 






Cross Section 




> 


<- 




Asymmetry - 


> 


INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





0. 


,2363 





.2364 


























1 


0. 


,0935 





.0935 


0. 


0935 


0. 


0935 






-0, 


.3521 


-0, 


.3521 


-0, 


.3521 


2 


0. 


,0954 





.0954 


0. 


0954 


0. 


0954 


0. 


,0954 


-0, 


.3522 


-0 


.3521 


-0, 


.3521 


3 


0. 


,0961 





.0961 


0. 


0961 


0. 


0961 






-0, 


.3525 


-0 


.3524 


-0, 


.3524 


4 


0. 


,4504 





.4504 


0. 


4505 










-0, 


.1802 


-0, 


.1802 






5 


0. 


,5480 





.5480 


0. 


5481 










-0, 


.0097 


-0 


.0097 






6 


0. 


,4414 





.4414 


0. 


4415 










-0, 


.1837 


-0, 


.1837 






7 


0. 


,5472 





.5472 


0. 


5473 










-0, 


.0097 


-0, 


.0097 






8 


0. 


,0000 





.0000 


0. 


0000 










, 


. 0000 


0, 


.0000 






9 


0. 


,5346 





.5346 


0. 


5347 










-0, 


.0103 


-0, 


.0103 






10 


2. 


,5217 


2 


.5217 










2. 


,5217 














11 


0. 


,4405 





.4404 


0. 


4404 


0. 


4406 






0, 


.6259 


0, 


.6258 


0, 


.6258 




<- 


— TauPol > 


























INDF 


ZUTPSM 


ZUTAU 


























3 


-0.0825 




-0.0826 
























SQRT(S) 


88.0000000 
















































<- 






Cross Section 




> 


<- 




Asymmetry - 


> 


INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





0. 


,3826 





.3827 


























1 


0. 


,1461 





.1461 


0. 


1461 


0. 


1461 






-0, 


.2723 


-0 


.2722 


-0, 


.2722 


2 


0. 


,1489 





.1489 


0. 


1490 


0. 


1490 


0. 


,1489 


-0, 


.2724 


-0 


.2723 


-0, 


.2723 


3 


0. 


,1500 





.1500 


0. 


1501 


0. 


1501 






-0, 


.2727 


-0 


.2726 


-0, 


.2726 


4 


0. 


,7092 





.7091 


0. 


7093 










-0, 


.1241 


-0, 


.1240 






5 


0. 


,8831 





.8831 


0. 


8832 










0, 


.0149 


0, 


.0149 






6 


0. 


,7000 





.6999 


0. 


7001 










-0, 


.1255 


-0 


.1255 






7 


0. 


,8824 





.8823 


0. 


8825 










0, 


.0149 


0, 


.0149 






8 


0. 


,0000 





.0000 


0. 


0000 










0, 


.0000 


0, 


.0000 






9 


0. 


,8635 





.8633 


0. 


8635 










0, 


.0145 


0, 


.0145 






10 


4. 


,0381 


4 


.0378 










4. 


,0381 














11 


0. 


,5124 





.5122 


0. 


5122 


0. 


5124 






0, 


.5569 


0, 


.5568 


0, 


.5569 



< TauPol > 

INDF ZUTPSM ZUTAU 
3 -0.0951 -0.0952 



73 



SQRKS) = 89.0000000 





<- 






Cross Section 




> 


<- 




Asymmetry - 


> 


INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





0. 


7007 


0, 


.7007 
























1 


0. 


2609 


0, 


.2608 


0. 


2609 


0.2609 






-0, 


,1870 


-0 


.1868 


-0, 


.1868 


2 


0. 


2658 


0, 


.2658 


0. 


2659 


0.2659 


0, 


.2658 


-0. 


.1870 


-0 


.1869 


-0, 


.1869 


3 


0. 


2677 


0, 


.2677 


0. 


2678 


. 2678 






-0. 


.1873 


-0 


.1871 


-0, 


.1871 


4 


1. 


2733 


1, 


.2731 


1. 


2734 








-0. 


,0662 


-0 


.0661 






5 


1. 


6126 


1, 


.6124 


1. 


6127 








0. 


,0388 





.0388 






6 


1. 


2634 


1, 


.2633 


1. 


2635 








-0. 


,0666 


-0 


.0666 






7 


1. 


6118 


1, 


.6117 


1. 


6120 










, uooo 





.0388 






8 


0. 


0000 


0, 


.0000 


0. 


0000 








0. 


,0000 





.0000 






9 


1. 


5792 


1, 


.5789 


1. 


5792 








0. 


,0385 





.0385 






10 


7. 


3403 


7, 


.3394 








7. 


.3402 














11 


0. 


6562 


0, 


.6558 


0. 


6561 


0.6562 






0. 


,4592 





.4591 


0, 


.4593 




<- 


■ — TauPol > 
























INDF 


ZUTPSM 


ZUTAU 
























3 


-0.1071 




-0.1072 






















SQRKS) 




90, 


.0000000 














































<- 






— Cross Section 




> 


<- 




Asymmetry - 


> 


INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





1. 


4868 


1, 


.4867 
























1 


0. 


5452 


0, 


.5451 


0. 


5453 


. 5453 






-0. 


.0993 


-0 


.0992 


-0, 


.0992 


2 


0. 


5555 


0, 


.5553 


0. 


5555 


0.5555 


0, 


.5554 


-0, 


.0994 


-0 


.0992 


-0, 


.0992 


3 


0. 


5593 


0, 


.5592 


0. 


5593 


0.5593 






-0. 


.0996 


-0 


.0994 


-0, 


.0994 


4 


2. 


6703 


2, 


.6698 


2. 


6704 








-0. 


.0086 


-0 


.0085 






5 


3. 


4173 


3, 


.4169 


3. 


4175 








0. 


,0617 





.0617 






6 


2. 


6586 


2, 


.6581 


2. 


6587 








-0. 


,0086 


-0 


.0085 






7 


3. 


4166 


3, 


.4161 


3. 


4167 








0. 


.0617 





.0617 






8 


0. 


0000 


0, 


.0000 


0. 


0000 








0, 


.0000 





.0000 






9 


3. 


3503 


3, 


.3495 


3. 


3501 








0. 


.0614 





.0615 






10 


15. 


5131 


15. 


.5104 








15. 


.5119 














11 


0. 


9678 


0, 


.9669 


0. 


9677 


0.9679 






0. 


.3302 





.3301 


0, 


.3303 



< TauPol > 

INDF ZUTPSM ZUTAU 
3 -0.1181 -0.1182 



74 



SQRT(S) = 91.0000000 



< Cross Section > < Asymmetry > 



INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





2, 


.8218 


2, 


.8214 
























1 


1, 


.0295 


1, 


.0292 


1 


.0295 


1, 


.0295 




-0. 


.0155 


-0 


.0153 


-0, 


.0153 


2 


1, 


.0484 


1 


.0481 


1 


.0485 


1, 


.0485 


1 . 0484 


-0, 


,0155 


-0 


.0153 


-0, 


.0153 


3 


1, 


.0556 


1, 


.0552 


1 


.0556 


1, 


.0556 




-0. 


.0156 


-0 


.0154 


-0, 


.0154 


4 


5, 


.0483 


5, 


.0472 


5 


.0483 








0. 


,0453 





.0454 






5 


6, 


.4869 


6 


.4858 


6 


.4870 








0. 


,0826 





.0827 






6 


5, 


.0335 


5, 


.0323 


5 


.0335 








0. 


,0454 





.0455 






7 


6, 


.4862 


6, 


.4850 


6 


.4862 








0. 


,0826 





.0827 






8 


0, 


.0000 





.0000 





.0000 








0. 


,0000 





.0000 






9 


6, 


.3635 


6, 


.3617 


6 


.3629 








0. 


,0825 





.0825 






10 


29, 


.4183 


29, 


.4120 










29.4150 














11 


1, 


.3601 


1 


.3581 


1 


.3600 


1, 


.3601 




0. 


.2053 





.2052 


0, 


.2054 




<- 


— TauPol > 
























INDF 


ZUTPSM 


ZUTAU 

























3 -0.1274 -0.1275 



SQRT(S) 
<- 



92.0000000 

Cross Section 



Asymmetry 



INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 





2. 


,3301 


2 


.3295 




















1 


0. 


,8526 





.8523 





.8525 


0. 


.8525 




0. 


.0530 





.0532 


2 


0. 


,8681 





.8677 





.8680 


0. 


.8680 


0.8681 


0. 


.0529 





.0531 


3 


0, 


,8739 





.8735 





.8738 


0. 


.8738 




0, 


,0529 





.0531 


4 


4. 


,1794 


4 


.1783 


4 


.1793 








0. 


,0884 





.0885 


5 


5. 


,3632 


5 


.3620 


5 


.3630 








0. 


,0993 





.0994 


6 


4. 


,1661 


4 


.1650 


4 


.1660 








0. 


,0886 





.0887 


7 


5. 


,3624 


5 


.3612 


5 


.3622 








0. 


,0994 





.0994 


8 


0. 


,0000 





.0000 





.0000 








0. 


,0000 





.0000 


9 


5. 


,2618 


5 


.2601 


5 


.2611 








0. 


,0993 





.0994 


10 


24, 


.3329 


24 


.3266 










24.3329 










11 


0. 


.9720 





.9701 





.9719 


0. 


.9719 




0. 


.1624 





.1624 




<- 


—TauPol > 




















INDF 


ZUTPSM 


ZUTAU 





















0.0532 
0.0531 
0.0531 



3 -0.1341 -0.1341 



75 



SQRT(S) = 93.0000000 



< Cross Section > < Asymmetry > 



INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





1. 


3333 


1. 


.3329 


























1 


0. 


4914 


0, 


,4912 





.4913 


0, 


,4913 






0, 


,1013 





.1015 


0. 


,1015 


2 


0. 


5002 


0, 


.5000 





.5001 


0. 


,5001 


0, 


,5002 


0. 


,1013 





.1015 


0. 


,1015 


3 


0. 


5035 


0. 


.5033 





.5035 


0. 


.5035 






0. 


,1013 





.1015 


0. 


,1015 


4 


2. 


4063 


2, 


.4056 


2 


.4062 










0. 


,1182 





.1183 






5 


3. 


0742 


3, 


.0734 


3 


.0740 










0. 


,1110 





.1111 






6 


2. 


3957 


2. 


.3950 


2 


.3956 










0. 


,1186 





.1187 






7 


3. 


0735 


3, 


.0727 


3 


.0733 










n 


1111 





.1111 






8 


0. 


0000 


0, 


.0000 





.0000 










0. 


,0000 





.0000 






9 


3. 


0156 


3. 


.0145 


3 


.0151 










0. 


,1111 





.1112 






10 


13. 


9653 


13, 


.9612 










13. 


.9658 














11 


0. 


5607 


0, 


.5595 





.5606 


0. 


.5606 






0. 


,2113 





.2113 


0. 


,2115 




<- 


■ — TauPol > 


























INDF 


ZUTPSM 


ZUTAU 


























3 


-0.1381 




-0.1380 
























SQRT(S) 




94, 


.0000000 
















































<- 






— Cross Section 




> 


<- 




Asymmetry - 


> 


INDF 


ZUTHSM 


ZUXSEC 


ZUXSA 


ZUXSA2 


ZUSMAT 


ZUTHSM 


ZUXSA 


ZUXSA2 





0. 


8343 


0, 


.8340 


























1 


0. 


3102 


0. 


.3100 





.3101 


0. 


.3101 






0. 


,1358 





.1360 


0. 


,1360 


2 


0. 


3157 


0, 


.3156 





.3157 


0, 


,3157 


0, 


.3157 


0, 


,1357 





.1359 


0. 


,1359 


3 


0. 


3178 


0, 


.3177 





.3178 


0. 


,3178 






0. 


,1357 





.1359 


0. 


,1359 


4 


1. 


5173 


1. 


.5168 


1 


.5172 










0. 


,1390 





.1391 






5 


1. 


9272 


1, 


.9267 


1 


.9270 










0. 


,1193 





.1194 






6 


1. 


5081 


1, 


.5076 


1 


.5080 










0. 


,1397 





.1398 






7 


1. 


9266 


1. 


.9260 


1 


.9264 










0. 


,1194 





.1194 






8 


0. 


0000 


0, 


.0000 





.0000 










0. 


,0000 





.0000 






9 


1. 


8899 


1, 


.8892 


1 


.8895 










0. 


,1195 





.1196 






10 


8. 


7692 


8. 


.7663 










8. 


.7692 














11 


0. 


3862 


0, 


.3854 





.3861 


0. 


,3862 






0. 


,3006 





.3006 


0. 


,3008 



< TauPol > 



INDF ZUTPSM ZUTAU 
3 -0 . 1404 -0 . 1403 



76 



SQRKS) = 95.0000000 

< Cross Section > 

INDF ZUTHSM ZUXSEC ZUXSA ZUXSA2 ZUSMAT 






0, 


.5808 


0. 


,5805 










1 


0, 


.2180 


0. 


,2178 


0, 


.2179 


0. 


.2179 


2 


0, 


.2218 


0. 


,2217 


0, 


.2218 


0. 


.2218 


3 


0, 


.2233 


0. 


,2232 


0, 


.2233 


0. 


.2233 


4 


1, 


.0652 


1. 


,0648 


1, 


.0650 






5 


1, 


.3442 


1. 


,3438 


1, 


.3441 






6 


1, 


.0567 


1. 


,0563 


1. 


.0566 






7 


1, 


.3436 


1. 


,3432 


1, 


.3434 






8 


0, 


.0000 


0. 


,0000 


0, 


.0000 






9 


1, 


.3177 


1. 


,3171 


1, 


.3174 






10 


6, 


.1274 


6. 


,1252 










11 


0, 


.3061 


0. 


,3056 


0, 


.3060 


0. 


.3061 




<- 


— TauPol > 











INDF ZUTPSM ZUTAU 
3 -0.1417 -0.1416 



< Asymmetry > 

ZUTHSM ZUXSA ZUXSA2 



0. 


.1617 


0. 


.1619 





0. 


.1615 


0. 


.1617 





0, 


,1615 


0, 


,1617 





0. 


,1542 


0. 


.1542 




0. 


,1255 


0. 


.1256 




0. 


,1552 


0, 


,1553 




0. 


,1256 


0. 


.1257 




0. 


.0000 


0. 


.0000 




0. 


,1258 


0. 


.1259 




0. 


.3902 


0. 


.3902 






.1619 
.1617 
.1617 



.3904 



77 



